Source: THE STANDARD HANDBOOK FOR AERONAUTICAL AND ASTRONAUTICAL ENGINEERS 


_ SECTION 3 _ 

MECHANICAL ENGINEERING 
PRINCIPLES 


Section Editor: Mark Davies 


3.1 

Downloaded from Digital Engineering Library @ McGraw-Hill (www.digitalengineeringlibrary.com) 
Copyright ©2004 The McGraw-Hill Companies. All rights reserved. 

Any use is subject to the Terms of Use as given at the website. 





MECHANICAL ENGINEERING PRINCIPLES 


PART 1 

STATICS OF RIGID BODIES 

Robert Paine 


In general, the study of mechanics may be divided into two distinct areas. These 
are statics, which involves the study of bodies at rest, and dynamics, which is the 
study of bodies in motion. In each case it is important to select an appropriate 
mathematical model from which a “free body diagram” may be drawn, representing 
the system in space, with all the relevant forces acting on that system. 

When a set of forces acts on a body, it gives rise to a resultant force or moment 

or a combination of both. The situation may be determined by considering three 
mutually perpendicular directions on the free body diagram and resolving the forces 
and moment in these directions. If the three directions are denoted by x, y, and z, 
then the sum of forces may be represented by 2F X , 2 F , and 2 F z and the sum of 
the moments about respective axes by 2 M x , 2 M , and 2 M z . Then for equilibrium 
the following conditions must hold: 

2F X = ZF y = J,F Z = 0 (3.1) 

2 M x = 2 M v = 2 M. = 0 (3.2) 

If the conditions in equations (3.1) and (3.2) are not satisfied, then there is a 
resultant force or moment, which is given by 

F = [(2 F x f + (2 F y f + (2 Ff] 112 

M = [(2M X ) 2 + (2M V ) 2 + (2M_) 2 ] 1/2 

The six conditions given in equations (3.1) and (3.2) satisfy problems in three 
dimensions. If one of these dimensions is not present (say, the z direction) the 
system reduces to a set of coplanar forces, and then 

2 F x = 2 M x = 2 M y = 0 

are automatically satisfied, and the necessary conditions of equilibrium in a two- 
dimensional system are 


2 F x = 2 F y = 2 M z = 0 (3.3) 

If the conditions in equation (1.3) are not satisfied, then the resultant force or 
moment is given by 

F = ((2 F x ) 2 + (2 F y ) 2 ] 1 ' 2 
M = 2 M. 

The above equations give solutions to what are said to be “statically determi¬ 
nate” systems. These are systems where there are the minimum number of con¬ 
straints to maintain equilibrium (Mathesun 1987). 


3.2 
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PART 2 

STRENGTH OF MATERIALS 

Robert Paine 


Weight: The weight ( W ) of a body is that force exerted due to gravitational attrac¬ 
tion on the mass ( m ) of the body: W = mg, where g is the acceleration due to 
gravity. 

Center of gravity: This is a point, which may or may not be within the body, 
at which the total weight of the body may be considered to act as a single force. 
The position of the center of gravity may be found experimentally or by analysis. 
When using analysis the moment of each element of weight, within the body, about 
a fixed axis is equated to the moment of the complete weight about that axis: 

x = HSing • x/1,Smg, y — 1,8>ng • yielding 

z — ?,Smg • z/'ZSmg 

where 8m is an element of mass at a distance of x, y or z from the respective axis, 
and x, y and z are the positions of the centers of gravity from these axes. Table 3.1 
shows the position of the center of gravity for some standard shapes. (See Meriam 
and Kraige 1987 for a more comprehensive list.) 

Shear force and bending moment: If a beam subject to loading, as shown in 
Figure 3.1, is cut, then in order to maintain equilibrium a shear force (Q) and a 
bending moment (M) must be applied to each portion of the beam. The magnitudes 
of Q and M vary with the type of loading and the position along the beam and are 
directly related to the stresses and deflections in the beam. 

Relationship between shearforce and bending moment: If an element of a beam 
is subjected to a load w then the following relationship holds: 

d 2 M d£ _ _ 
dx 2 dx W 

Table 3.2 shows examples of bending moments, shear force, and maximum deflec¬ 
tion for standard beams. 

Bending equation: If a beam has two axes of symmetry in the xy plane then the 
following equation holds: 


Mzfh ~ EIR z = cr/y 

where M z is the bending moment, R. the radius of curvature, L the moment of 
inertia, E the modulus of elasticity, y the distance from the principal axis, and cr 
the stress. 

Torsion equation: If a circular shaft is subject to a torque ( T ) then the following 
equation holds: 


777 = rlr = G6/L 

where 7 is the polar second moment of area, G the shear modulus, L the length, 0 
the angle of twist, t the shear stress, and r the radius of the shaft. 


3.3 
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3.4 SECTION THREE 

TABLE 3.1 Centers of Gravity and Moments of Inertia or Second Moments of Area for 
Two-Dimensional Figures 

Shape G I 

y = hl 3 I GG = bh 3 /36 

l xx = bh 3 /n 


Triangular area 


Rectangular-area 



Circular sector 



x 


2r sin a 
3 a 


Slender rod 


Rectangular prism 


Thin disk 


l 




Circular cylinder 



I GG = Wi 3 /12 
l xx = bhV 3 



Aw ^77 ml 2 / 12 


I xx = m(b 2 + c 2 )/ 12 
I Y y = m(c 2 + a 2 )/12 
I zz = m(a 2 + b 2 )/12 

/ XA; = mr 2 /2 
try = 4z = '»f 2 A 4 


Atjt = f>ir 2 !2 

lyy = 4z = m(3r 2 + L 2 )I12 


Sphere 


1 w fry fzz 2 tnr 2 / 5 


Hemisphere 

Circular cone 



y = 3r/8 

x = h / 4 


lyy = 2mr 2 /5 

I xx = 3nir 2 /10 

3mr 2 m/t 2 

fyy — - "f - 

" 20 10 
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FIGURE 3.1 
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3.6 SECTION THREE 

TABLE 3.2 



w 


A 


B 

f 


mo 



'Wi 2 

2 IT 


w 


|C A 


B D[ 

1 


L 

2 

2 

¥ 



One concentrated load W 
M at A = Wx, Q at A = W 
M greatest at B, and = WL 
Q uniform throughout 

Maximum deflection = WL 3 /3EI at the free end. 


Uniform load of W 
M at A = Wx 2 /2L 
Q at A = Wx/L 
M greatest at B = WL/2 
Q greatest at B = W 

Maximum deflection = WL2/8EI at the free end. 
One concentrated load at the center of a beam 


W 

M at A = — 
2 


Q at A = W/2 

M greatest at B = WL/A 
Q uniform throughout 

Maximum deflection = WL 3 /48 El at the center 
Uniform load W 


W L 2 


M at A = —-x 

2L V 4 


Q at A = Wx/L 
M greatest at B = WL/ 8 

Q greatest at C and D = W/2 maximum deflection at 
B = 5WX 3 /384£7 

Beam fixed at ends and loaded at center. 

M is maximum at A , B, and C and = WL/S. 
Maximum deflection at C = WX 3 /192£7 


Beam fixed at ends with uniform load. 
M maximum at A and B and = WL /12 
Maximum deflection at C = WX 3 /384£7 


One concentrated load W 
Reaction R = 5VK/16 
M maximum at A, and = 3WX/16 
M at C = 5WL/32 

Maximum deflection is L/V5 from the free end, and 

= WL 3 /101EI 

Uniform load W 

Reaction R = 3VV78 

M maximum at A , and = WL/ 8 

M at C = 9VKL/128 

Maximum deflection is 3Z./8 from the free end, and = 
WL 2 /1Z1EI 
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MECHANICAL ENGINEERING PRINCIPLES 3.7 

PART 3 

DYNAMICS OF RIGID BODIES 

Robert Paine 


3.1 BASIC DEFINITIONS 


Newton's Laws of Motion 

First Law. A particle remains at rest or continues to move in a straight line with 
a constant velocity unless acted on by an external force. 

Second Law. The sum of all the external forces acting on a particle is propor¬ 
tional to the rate of change of momentum. 

Third Law. The forces of action and reaction between interacting bodies are 
equal in magnitude and opposite in direction. 

Newton’s law of gravitation, which governs the mutual interaction between bod¬ 
ies, states: 


F = Gm l m 2 /x 2 

where F is the mutual force of attraction, G is a universal constant called the 
constant of gravitation, which has a value 6.673 X 10~ n m 3 kg -1 s 2 , m x and m 2 are 
the masses of the two bodies, and x is the distance between the centers of the 
bodies. 

Mass (m) is a measure of the amount of matter present in a body. 

Velocity is the rate of change of distance (x) with time (t): 

v = dxl dt or x 

Acceleration is the rate of change of velocity ( v ) with time ( t ): 

a = dv I dt or d 2 x/dt 2 or x 

Momentum is the product of the mass and the velocity. If no external forces are 
present then the momentum of any system remains constant. This is known as the 
conservation of momentum. 

Force is equal to the rate of change of momentum (mv) with time (t): 

F = d{mv)/dt 
F = m ■ dv/dt + v ■ dmldt 

If the mass remains constant then this simplifies to F = m ■ dv/dt, i.e., Force = 
mass X acceleration, and it is measured in newtons. 

Impulse (/) is the product of the force and the time that force acts. Since I = 
Ft = mat = m(v 2 — iq), impulse is also said to be the change in momentum. 

Energy: There are several different forms of energy which may exist in a system. 
These may be converted from one type to another but they can never be destroyed. 
Energy is measured in joules. 
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3.8 SECTION THREE 


Potential energy (PE) is the energy which a body possesses by virtue of its 
position in relation to other bodies: PE = mgh, where h is the distance above some 
fixed datum and g is the acceleration due to gravity. 

Kinetic energy (KE) is the energy a body possesses by virtue of its motion: 
KE = Vimv 2 . 

Work (W) is a measure of the amount of energy produced when a force moves 
a body a given distance: W — F • x. 

Power (P ) is the rate of doing work with respect to time and is measured in 
watts. 

Moment of inertia (I): The moment of inertia is that property in a rotational 
system which may be considered equivalent to the mass in a translational system. 
It is defined about an axis xx as l xx = Y,8mx 2 = mk 2 XX , where x is the perpendicular 
distance of an element of mass 8m from the axis xx and k xx is the radius of gyration 
about the axis xx. Table 3.1 gives some data on moments of inertia for standard 
shapes. 

Angular velocity (co) is the rate of change of angular distance (6) with time: 

= dO/dt — 8 

Angular acceleration (cr) is the rate of change of angular velocity (co) with time: 


= dco/dt or d 2 9/dt 2 or 6 

Both angular velocity and accleration are related to linear motion by the equations 
v = cox and a = ocx (see Figure 3.2). 

Torque (T) is the moment of force about the axis of rotation: 

T — IqOc 
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A torque may also be equal to a couple, which is two forces equal in magnitude 
acting some distance apart in opposite directions. 

Parallel axis theorem: if I GG is the moment of inertia of a body of mass m about 
its center of gravity, then the moment of inertia (I) about some other axis parallel 
to the original axis is given by I = I GG + mr 2 , where r is the perpendicular distance 
between the parallel axes. 

Perpendicular axis theorem: If l xx , l Yy , and l /:z represent the moments of inertia 
about three mutually perpendicular axes x, y, and z for a plane figure in the xy 
plane (see Figure 3.3) then I /z — l xx + I ry . 

Angular momentum ( H a ) of a body about a point O is the moment of the linear 
momentum about that point and is col 00 . The angular momentum of a system re¬ 
mains constant unless acted on by an external torque. 

Angular impulse is the product of torque by time, i.e. angular impulse = Tt = 
la ■ t = I(co 2 — &>,), the change in angular momentum. 

Angular kinetic energy about an axis O is given by ViIqCo 2 . 

Work done due to a torque is the product of torque by angular distance and is 
given by T6. 

Power due to torque is the rate of angular work with respect to time and is 
given by TdO!dt = Too. 

Friction: Whenever two surfaces, which remain in contact, move one relative to 
the other there is a force which acts tangentially to the surfaces so as to oppose 
motion. This is known as the force of friction. The magnitude of this force is pR, 
where R is the normal reaction and p is a constant known as the coefficient of 
friction. The coefficient of friction depends on the nature of the surfaces in contact. 


3.2 LINEAR AND ANGULAR MOTION IN TWO 
DIMENSIONS 


Constant acceleration: If the accleration is integrated twice and the relevant initial 
conditions are used, then the following equations hold: 
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3.10 


SECTION THREE 


Linear motion Angular motion 

x = v{t + Viat 2 9 = cop + Viat 2 

V2 = Vi Clt (x>2 = Ct>i + CXt 

v\ ~ v i + 2ax a> 2 = co j + 2 ad 


Variable acceleration: If the acceleration is a function of time then the area 
under the acceleration time curve represents the change in velocity. If the acceler¬ 
ation is a function of displacement then the area under the acceleration distance 
curve represents half the difference of the square of the velocities (see Figure 3.4). 
Curvilinear motion is when both linear and angular motions are present. 

If a particle has a velocity v and an acceleration a then its motion may be 
described in the following ways: 

1. Cartesian components, which represent the velocity and acceleration along two 
mutually perpendicular axes x and y (see Figure 3.5(a)): 


v x = v cos 9, v Y — v sin 6 , a x = a cos <f> 


Cly 


a sin 4 > 


a 


a 


2 



t 


x 


a = —- or atit - dv 
df 


a 


dx . dv 
df dx 


Area 


r adf 


dt’ 

dx 


= V 2 -Vi 


or adx = vdv 



FIGURE 3.4 
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FIGURE 3.5 


2. Normal and tangential components: see Figure 3.5(b): 

v, = v = rb = rco, v„ — 0 
a, = rb + rb = ra + rco 
a n = vb = rco 2 

3. Polar coordinates: see Figure 3.5(c): 

v, — r, v e — rb 

a, = r — rb 2 , a e = rb + 2 rb 

3.3 CIRCULAR MOTION 


Circular motion is a special case of curvilinear motion in which the radius of 
rotation remains constant. In this case there is an acceleration towards the center 
of co 2 r. This gives rise to a force towards the center known as the centripetal force. 
This force is reacted to by what is called the centrifugal reaction. 

Velocity and acceleration in mechanisms: A simple approach to determine the 
velocity and acceleration of a mechanism at a point in time is to draw velocity and 
acceleration vector diagrams. 

Velocities: If in a rigid link AB of length / the end A is moving with a different 
velocity to the end B, then the velocity of A relative to B is in a direction perpen¬ 
dicular to AB (see Figure 3.6). 

When a block slides on a rotating link the velocity is made up of two compo¬ 
nents, one being the velocity of the block relative to the link and the other the 
velocity of the link. 

Accelerations: If the link has an angular acceleration a then there will be two 
components of acceleration in the diagram, a tangential component al and a cen¬ 
tripetal component of magnitude co 2 l acting towards A. 

When a block slides on a rotating link the total acceleration is composed of four 
parts: first, the centripetal acceleration towards O of magnitude orl; second, the 
tangential acceleration ah, third, the acceleration of the block relative to the link; 
fourth, a tangential acceleration of magnitude 2 vco known as Coriolis acceleration. 
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3.12 


SECTION THREE 




ed 

(b) 

FIGURE 3.6 


The direction of Coriolis acceleration is determined by rotating the sliding velocity 
vector through 90° in the direction of the link angular velocity o>. 


3.4 LINEAR AND ANGULAR MOTION IN THREE 
DIMENSIONS 


Motion of a Particle in a Moving Coordinate System 

xyz is a moving coordinate system, with its origin at O, which has a position vector 
R. a translational velocity vector R and an angular velocity vector oj relative to a 
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fixed coordinate system XYZ, origin at O'. Then the motion of a point P whose 
position vector relative to O is P and relative to O' is r is given by the following 
equations (see Figure 3.7): 


r = R + p r + toXp 

where p,. is the velocity of the point P relative to the moving system xyz and to X 
p is the vector product of to and p: 

r = R + toXp + toX(toXp) + 2to x p, + p, 

where p r is the acceleration of the point P relative to the moving system. Thus r 
is the sum of: 

1. The relative velocity p ( . 

2. The absolute velocity R of the moving origin O 

3. The velocity to X p due to the angular velocity of the moving axes xyz 
and r is the sum of: 

1. The relative acceleration p r 

2. The absolute acceleration R of the moving origin O 

3. The tangential acceleration to X p due to the angular acceleration of the moving 
axes xyz 

4. The centripetal acceleration to X (to X p) due to the angular velocity of the 
moving axes xyz 

5. Coriolis component acceleration 2to X p r due to the interaction of coordinate 
angular velocity and relative velocity 
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In all the vector notation a right-handed set of coordinate axes and the right-hand 
screw rule is used. 


Gyroscopic Effects 

Consider a rotor which spins about its geometric axis (see Figure 3.8) with an 
angular velocity to. Then two forces F acting on the axle to form a torque T, whose 
vector is along the x-axis, will produce a rotation about the y axis. This is known 
as precession, and it has an angular velocity Cl. It is also the case that if the rotor 
is precessed then a torque T will be produced, where T is given by T = l xx (oil. 
When this is observed it is the effect of gyroscopic reaction torque that is seen, 
which is in the opposite direction to the gyroscopic torque (Scarborough 1958). 


3.5 BALANCING 


In any rotational or reciprocating machine where accelerations are present, unbal¬ 
anced forces can lead to high stresses and vibrations. The principle of balancing is 
such that by the addition of extra masses to the system the out-of-balance forces 
may be reduced or eliminated. 


3.6 BALANCING OF ROTATING MASSES 


Single Out-of-Balance Mass 

One mass (m) at a distance r from the center of rotation and rotating at a constant 
angular velocity to produces a force mto 2 r. This can be balanced by a mass M 
placed diametrically opposite at a distance R, such that MR = mr. 


V 



Precession axis 
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Several Out-of-Balance Masses in One Transverse Plane 

If a number of masses (m u m 2 , . . .) are at radii (ry, r 2 , . . .) and angles (8 lt 0 2 , 
. . .) (see Figure 3.9), then the balancing mass M must be placed at a radius R such 
that MR is the vector sum of all the mr terms. 


Masses in Different Transverse Planes 

If the balancing mass in the case of a single out-of-balance mass were placed in a 
different plane then the centrifugal force would be balanced. This is known as static 
balancing. However, the moment of the balancing mass about the original plane 
would lead to what is known as dynamic unbalance. 

To overcome this, the vector sum of all the moments about the reference plane 
must also be zero. In general, this requires two masses placed in convenient planes 
(see Figure 3.10). 


Balancing of Reciprocating Masses in Single-Cylinder Machines 

The acceleration of a piston as shown in Figure 3.11 may be represented by the 
equation (Wilson 1959) 

x = —&> 2 r[cos 8 + (l/«)cos 28 + (l/4n) (cos 28 — cos 49) + ....* 

where n — Mr. If n is large, then the equation may be simplified and the force 
given by 



£A X = Xmto 2 r sin 6 = 0 
2 F y = cos 0 = 0 

FIGURE 3.9 


- This equation forms an infinite series in which higher terms are small, and they may be ignored for 
practical situations. 
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X Ax = Xmu 2 r sin 6 = 0 and XFy = Xmco 2 r cos 0 = 0 
as in the previous case, also 
XM x = XmcFr sin 6 . a = 0 
X My = Xmcfr cos 6 . a = 0 

FIGURE 3.10 



W M 

FIGURE 3.11 


F = mx — —m<u 2 r[cos 6 + (l/n)cos 26] 

The term mco 2 r cos 6 is known as the primary force and (1 /n)ma> 2 r cos 26 as the 
secondary force. Partial primary balance is achieved in a single-cylinder machine 
by an extra mass M at a radius R rotating at the crankshaft speed. Partial secondary 
balance could be achieved by a mass rotating at 2a>. As this is not practical this is 
not attempted. When partial primary balance is attempted a transverse component 
Mit> 2 R sin 6 is introduced. The values of M and R are chosen to produce a com¬ 
promise between the reciprocating and the transverse components. 


Balancing of Reciprocating Masses in Multicylinder Machines 

When considering multicylinder machines, account must be taken of the force pro¬ 
duced by each cylinder and the moment of that force about some datum. The 
conditions for primary balance are 

F = 'Zmorir cos 0 = 0, M = 1,mco 2 r cos 6 ■ a = 0 

where a is the distance of the reciprocating mass m from the datum plane. 
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In general, the cranks in multicylinder engines are arranged to assist primary 
balance. If primary balance is not complete, then extra masses may be added to the 
crankshaft but these will introduce an unbalanced transverse component. The con¬ 
ditions for secondary balance are 

F = 'Emio 2 (r/n) cos 20 = Sw(2tu) 2 (r/4«) cos 20 = 0 


and 


M = Y.m(2w) 2 (r/4n) cos 20 ■ a — 0 

The addition of extra masses to give secondary balance is not attempted in practical 
situations. 
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PART 4 

VIBRATIONS 

Robert Paine 
Christopher Beards 


3.7 SINGLE-DEGREE-OF-FREEDOM SYSTEMS 


The term degrees of freedom in an elastic vibrating system is the number of pa¬ 
rameters required to define the configuration of the system. To analyze a vibrating 
system, a mathematical model is constructed, which consists of springs and masses 
for linear vibrations. The type of analysis then used depends on the complexity of 
the model. 

Rayleigh’s method: Rayleigh showed that if a reasonable deflection curve is 
assumed for a vibrating system, then by considering the kinetic and potential en¬ 
ergies* an estimate to the first natural frequency could be found. If an inaccurate 
curve is used then the system is subject to constraints to vibrate it in this unreal 
form, and this implies extra stiffness such that the natural frequency found will 
always be high. If the exact deflection curve is used then the natural frequency will 
be exact. 


Transverse Vibration of Beams 

Consider a beam of length (/), weight per unit length (w), modulus (£), and moment 
of inertia (/). Then its equation of motion is given by 


* Consider the equation of motion for an undamped system (Figure 3.13): 


d 2 x 

m • —- + kx = 0 
dt 2 


(3.4) 


but 


d 2 x d (dx\ dx d (dx\ 1 d (dx 

dt 2 dt \dt) dt dx\dt) 2 dx\dt 


Therefore, equation (3.4) becomes 


d 

dx 


1 / dx 
- m\ — 

2 \dt 


+ kx = 0 


Integrating gives 


1 dxY 1 0 

- m — H— kx 2 = Constant 

2 \dt 2 


the term Vim(dxl dt) 2 represents the kinetic energy and Vikx 2 the potential energy. 


3.18 
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d^y 
dx 4 


El-fr* ~ W(o 2 y/g = 0 


where co is the natural frequency. The general solution of this equation is given by 
y = A cos fix + B sin fix + C cos h fix + D sin h fix 
where ff = wco 2 /gEI. 

The four constants of integration A, B, C, and D are determined by four inde¬ 
pendent end conditions. In the solution trigonometrical identities are formed in /3 
which may be solved graphically, and each solution corresponds to a natural fre¬ 
quency of vibration. Table 3.3 shows the solutions and frequencies for standard 
beams (Young and Feglar 1949). 

Dunkerley’s empirical method is used for beams with multiple loads. In this 
method the natural frequency (/[) is found due to just one of the loads, the rest 
being ignored. This is repeated for each load in turn and then the natural frequency 
of vibration of the beam due to its weight alone is found (/ 0 ). Then the natural 
frequency of vibration of the complete system (/) is given by 


J_ J_ J_ J_ J_ 

f2 ~ f2 -C2 -T2 f2 ~*~ 

J Jo J 1 J 2 J 3 


1 




2 

n 


(see Cole 1950 for a more detailed explanation). 

Whirling of shafts: If the speed of a shaft or rotor is slowly increased from rest 
there will be a speed where the deflection increases suddenly. This phenomenon is 
known as whirling. Consider a shaft with a rotor of mass m such that the center of 
gravity is eccentric by an amount e. If the shaft now rotates at an angular velocity 
co then the shaft will deflect by an amount y due to the centrifugal reaction (see 
Figure 3.12). Then 


mco 2 (y + e) = ky 

where k is the stiffness of the shaft. Therefore, 

_ e 

(I/mo 2 — 1) 

When ( k/ma> 2 ) = 1, y is then infinite and the shaft is said to be at its critical 
whirling speed co c . At any other angular velocity co the deflection y is given by 



When co < w c , y is the same sign as e and as to increases towards co c the deflection 
theoretically approaches infinity. When co > co c , y is opposite in sign to e and will 
eventually tend to —e. This is a desirable running condition with the center of 
gravity of the rotor mass on the static deflection curve. Care must be taken not to 
increase co too high, as co might start to approach one of the higher modes of 
vibration (Thompson 1983). 

Torsional vibrations: The following section deals with transverse vibrating sys¬ 
tems with displacements x and masses m. The same equations may be used for 
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torsional vibrating systems by replacing x by 9 the angular displacement and m by 
/, the moment of inertia. 


Undamped Free Vibrations 

The equation of motion is given by mx + kx = 0 or jc + oj; r x = 0, where m is the 
mass, k the stiffness and cof, = klm, which is the natural frequency of vibration of 
the system (see Figure 3.13). The solution to this equation is given by 

x = A sin(cu„f + a) 

where A and a are constants which depend on the initial conditions. This motion 
is said to be simple harmonic with a time period T = 2irlu> n . 


Damped Free Vibrations 

The equation of motion is given by mx + cx + kx = 0 (see Figure 3.14), where c 
is the viscous damping coefficient, or Jc + ( chn)x + co^x = 0. The solution to this 
equation and the resulting motion depends on the amount of damping. If c > 2mco n 
the system is said to be overdamped. It will respond to a disturbance by slowly 
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FIGURE 3.13 


FIGURE 3.14 


'77777777/ 

X 


returning to its equilibrium position. The time taken to return to this position de¬ 
pends on the degree of damping (see Figure 3.15(c)). If c = 2 mco n the system is 
said to be critically damped. In this case it will respond to a disturbance by returning 
to its equilibrium position in the shortest possible time. In this case (see Figure 
3.15(b)) 


—— g—(cl2m)t(A+Bt) 

where A and B are constants. If c < 2 ma> n , the system has a transient oscillatory 
motion given by 

x = e -( c/2m > r [C sin(tu^ - c 2 /4m 2 ) U2 t + D( cos a> 2 - c 2 /4m 2 ) 1/2 t\ 

where C and D are constants. The period 

T 2t7 

= (co 2 - c 2 /4m 2 ) 1/2 


(see Figure 3.15(a)). 

Logarithmic Decrement. A way to determine the amount of damping in a system 
is to measure the rate of decay of successive oscillations. This is expressed by a 
term called the logarithmic decrement (S), which is defined as the natural logarithm 
of the ratio of any two successive amplitudes (see Figure 3.16): 

S = log r .(x,/x 2 ) 


where x is given by 




+ <l> 


Therefore, 
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S = log (c < ' < / ' 

= cr/2m 

where r is the period of damped oscillation. 

If the amount of damping present is small compared to the critical damping, r 
approximates to 2 v! co, and then 


S = CTr/mco n 


Forced Undamped Vibrations 

The equation of motion is given by (see Figure 3.17) 

mx + kx = F 0 sin cot 


or 


x + of; = ( F 0 /m ) sin cot 

The solution to this equation is 

x = C sin co n t + D cos co n t + F 0 cos cot\m(col — to 2 )] 

where co is the frequency of the forced vibration. The first two terms of the solution 
are the transient terms, which die out, leaving an oscillation at the forcing frequency 
of amplitude 
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3.25 


F 0 /\m(aol - to 2 )] 


or 


Fp ( <*>l 

k — co 2 

The term toll {to 2 — co 2 ) is known as the dynamic magnifier, and it gives the 
ratio of the amplitude of the vibration to the static deflection under the load F 0 . 
When to = to n the amplitude becomes infinite and resonance is said to occur. 


Forced Damped Vibrations 

The equation of motion is given by (see Figure 3.17(b)) 

mx + cx + kx = F 0 sin cot 


or 


x + (clm)x + co 2 = ( F 0 /m ) sin cot 

The solution to this equation is in two parts: a transient part as in the undamped 
case which dies away, leaving a sustained vibration at the forcing frequency given 
by 


X 


Fc _ 

m [(to 2 


1 

co 2 ) 2 + (cdWm) 2 ] 1 ' 2 


sin(cuf 


a ) 


The term 


[{co 2 n ~ " 2 ) 2 + (cWm)] 2 ] 1/2 

is called the dynamic magnifier. Resonance occurs when co — co n . As the damping 
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is increased the value of cu for which resonance occurs is reduced. There is also a 
phase shift as w increases tending to a maximum of ir radians. It can be seen in 
Figure 3.18(a) that when the forcing frequency is high compared to the natural 
frequency the amplitude of vibration is minimized. 


Forced Damped Vibrations Due to Reciprocating or Rotating Unbalance 

Figure 3.19 shows two elastically mounted systems, (a) with the excitation supplied 
by the reciprocating motion of a piston, and (b) by the rotation of an unbalanced 
rotor. In each case the equation of motion is given by 




Frequency ratio (cj/co n ) 
(b) 


FIGURE 3.18 
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(M — m)x + cx + Lx = (mew 2 ) sin tot 


The solution of this equation is a sinusoid whose amplitude, X, is given by 


meto 2 


X = 


V[(^ - Mco 2 ) 2 + (cut) 2 ] 


In representing this information graphically it is convenient to plot MX /me against 
to/co n for various levels of damping (see Figure 3.20(a)). From this figure it can be 
seen that for small values of to the displacement is small, and as to is increased the 
displacement reaches a maximum when to is slightly greater than to n . As to is further 
increased the displacement tends to a constant value such that the center of gravity 
of the total mass M remains stationary. Figure 3.20(b) shows how the phase angle 
varies with frequency. 

Forced Damped Vibration. If a system as shown in Figure 3.21 has a sinusoidal 
displacement applied to its base of amplitude, y, then the equation of motion be¬ 
comes 


mx + cx + kx = ky + cy 


The solution of this equation yields 



where x is the amplitude of motion of the system. 

When this information is plotted as in Figure 3.22, it can be seen that for very 
small values of to the output amplitude X is equal to the input amplitude Y. As to 
is increased towards to,, the output reaches a maximum. When to = V2 to,, the 
curves intersect and the effect of damping is reversed. 

The curves in Figure 3.22 may also be used to determine the amount of sinu¬ 
soidal force transmitted through the springs and dampers to the supports, i.e., the 
axis (X/Y) may be replaced by (F r /F 0 ) where F 0 is the amplitude of applied force 
and F t is the amplitude of force transmitted. 


m 


x 


'//////, 


c 


k 


V 


'////////////////////A 


FIGURE 3.21 
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(a) 



FIGURE 3.22 


Frequency ratio (oj/cc n ) 
(b) 


3.8 MULTI-DEGREE-OF-FREEDOM SYSTEMS 


Normal Mode Vibration 

The fundamental techniques used in modeling multi-degree-of-freedom systems 
may be demonstrated by considering a simple two-degree-of-freedom system as 
shown in Figure 3.23. The equations of motion for this system are given by 


m 1 .x l + (k 1 + k 2 )x 1 — k 2 x 2 = 0 


m 2 x 2 + (k 2 + k 2 )x 2 — k 2 x 1 = 0 


or in matrix form: 
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Assuming the motion of every point in the system to be harmonic, then the 
solutions will take the form 


x 1 = sin cot 
x 2 = A 2 sin cot 

where A l and A 2 are the amplitudes of the respective displacements. By substituting 
the values of x v x 2 , x v and x 2 into the original equations the values of the natural 
frequencies of vibration may be found along with the appropriate mode shapes. 
This is a slow and tedious process, especially for systems with large numbers of 
degrees of freedom, and is best performed by a computer program. 


The Holtzer Method 

When only one degree of freedom is associated with each mass in a multimass 
system, then a solution can be found by proceeding numerically from one end of 
the system to the other. If the system is being forced to vibrate at a particular 
frequency, then there must be a specific external force to produce this situation. A 
frequency and a unit deflection is assumed at the first mass and from this the inertia 
and spring forces are calculated at the second mass. This process is repeated until 
the force at the final mass is found. If this force is zero, then the assumed frequency 
is a natural frequency. Computer analysis is most suitable for solving problems of 
this type. 

Consider several springs and masses as shown in Figure 3.24. Then with a 
unit deflection at the mass m 1 and an assumed frequency co there will be an inertia 



FIGURE 3.24 
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force of w,a> 2 acting on the spring with stiffness k 1 . This causes a deflection of 
m x (i> 2 lk x , but if m 2 has moved a distance x 2 then m l (xi 2 lk l = 1 — x 2 or x 2 = 
1 — The inertia force acting due to m 2 is m 2 (o 2 x 2 , thus giving the total 

force acting on the spring of stiffness k 2 as [mj® 2 + m 2 co 2 x 2 \/k 2 . Hence the dis¬ 
placement at x 3 can be found and the procedure repeated. The external force acting 
on the final mass is then given by 


n 

X wt,.o> 2 x,.. 
1 = 1 


If this force is zero, then the assumed frequency is a natural one. 


3.9 RANDOM VIBRATIONS 


Introduction 

If the vibration response parameters of a dynamic system are accurately known as 
functions of time, the vibration is said to be deterministic. However, in many sys¬ 
tems and processes responses cannot be accurately predicted; these are called ran¬ 
dom processes. Examples of a random process are turbulence, fatigue, the meshing 
of imperfect gears, surface irregularities, the motion of a car running along a rough 
road, and building vibration excited by an earthquake (Figure 3.25). 

A collection of sample functions x t (t), x 2 (t), x 3 (t), . . . , x n (f) which make up the 
random process x(t ) is called an ensemble (Figure 3.26). These functions may com¬ 
prise, for example, records of pressure fluctuations or vibration levels, taken under 
the same conditions but at different times. 

Any quantity which cannot be precisely predicted is nondeterministic and is 
known as a random variable or a probabilistic quantity. That is, if a series of tests 
are conducted to find the value of a particular parameter, x, and that value is found 
to vary in an unpredictable way that is not a function of any other parameter, then 
x is a random variable. 



FIGURE 3.25 Example random process variable as /(f). 
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* n (r) 





FIGURE 3.26 Ensemble of a random process. 


Probability Distribution 

If n experimental values of a variable x are x 1; x 2 , x 3 , . . . . x n , the probability that 
the value of x will be less than x' is n'In, where n' is the number of x values which 
are less than or equal to x'. That is, 

Prob(.r =s x') = n' In 

When n approaches °o this expression is the probability distribution function of x, 
denoted by P(x), so that 

P{x) = Lt (n'In) 

n —►co 

The typical variation of P{x) with x is shown in Figure 3.27. Since x(t) denotes a 
physical quantity, 



FIGURE 3.27 Probability distribution function as fix). 
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Prob(x < — go) = 0, and Prob(x < + oo) = 1 

The probability density function is the derivative of P(x) with respect to x and 
this is denoted by p(x). That is, 


p(x) 


dPjx) 

dx 


= Lt 

Ax—*0 


P(x + Ax) 
Ax 


P(x) 


where P(x + Ax) — Pix) is the probability that the value of x(t) will lie between x 
and x + Ax (Figure 3.27). Now 


so that 


Hence, 


p(x) 


dP(x) 

dx 


Pix) 



p{x)dx 


P(cc) = p(x)dx — 1 


so that the area under the probability density function curve is unity. 

A random process is stationary if the joint probability density 

Pix){t i), x(t 2 ), x(t 3 ), . . .) 

depends only on the time differences t 2 — t x , t 3 — t 2 , and so on, and not on the 
actual time instants. That is, the ensemble will look just the same if the time origin 
is changed. A random process is ergodic if every sample function is typical of the 

entire group. _ 

The expected value of /(x), which is written E[f(x)] or f(x) is 

£T.f(x)l = fix) = f(x)p(x)dx 

so that the expected value of a stationary random process x(t) is 

E[xif)] = E[x(t t + t )] 


for any value of t. 

If f{x) = x, the expected value or mean value of x, £[x] or x, is 
E[x\ = x = J xp(x)dx 

In addition, if fix) = x 2 , the mean square value of x, x 2 is 
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E[x 2 ~\ = x 2 — J x 2 ■ p(x)dx 


The variance of x, cr 2 is the mean square value of x about the mean, that is. 



cr is the standard deviation of x, hence 


Variance = (Standard deviation) 2 

= {Mean square — (Mean) 2 } 


If two or more random variables x, and x 2 represent a random process at two 
different instants of time, then 



and if t 2 and t 2 are the two instants of time, 

E[x(t,),x(t 2 )\ = R(t,,t 2 ) 

which is the autocorrelation function for the random process (Figure 3.28). 

For random processes which are stationary, 

E[x(t\),x(t 2 j\ = R(t\,t 2 ) — E(t 2 — h) = R(t), say, 

since the average depends only on time differences. If the process is also ergodic, 
then 


Xnlt) 


X 3 (t) 


*l( 



FIGURE 3.28 Random processes. 
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It is worth noting that 


r( o) = Eum 



x 2 (t)dt 


which is the average power in a sample function. 


Random Processes 

The most important random process is the Gaussian or normal random process. 
This is because a wide range of physically observed random waveforms can be 
represented as Gaussian processes, and the process has mathematical features which 
make analysis relatively straightforward. 

The probability density function of a Gaussian process x{t) is 



where cr is the standard deviation of x and x is the mean value of x. The values of 
cr and x may vary with time for a non-stationary process but are independent of 
time if the process is stationary. 

One of the most important features of the Gaussian process is that the response 
of a linear system to this form of excitation is usually another (but still Gaussian) 
random process. The only changes are that the magnitude and standard deviation 
of the response may differ from those of the excitation. 

A Gaussian probability density function is shown in Figure 3.29. It can be seen 
to be symmetric about the mean value x, and the standard deviation cr controls the 
spread. 

The probability that x{t) lies between — Acr and +Acr, where A is a positive 
number, can be found since, if x = 0, 


1 


— i — e Mii^i„h dx 
A<r \Z2tt cr 


Prob{-Acr — x(t) -- + Ao-} 


Figure 3.30 shows the Gaussian probability density function with zero mean. This 
integral has been calculated for a range of values of A and the results are given in 
Table 3.4. The probability that x(t) lies outside the range — Ao- to +Ao- is 1 minus 
the value of the above integral. This probability is also given in Table 3.4. 


Spectral Density 

The spectral density S(co) of a stationary random process is the Fourier transform 
of the autocorrelation function R(t), and is given by 



The inverse, which also holds true, is 
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FIGURE 3.29 Gaussian probability density function. 



FIGURE 3.30 Gaussian probability density function with zero 
mean. 


R(t) = J S(co)e ,01T dco 

If r = 0 

R( 0) = f S((o)dw 
= E[x 2 ] 

That is, the mean square value of a stationary random process x is the area under 
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TABLE 3.4 


Value of A 

Prob[— Act - x(t ) *£ Act] 

Prob[|x(f)| > A or] 

0 

0 

1.0000 

0.2 

0.1585 

0.8415 

0.4 

0.3108 

0.6892 

0.6 

0.4515 

0.5485 

0.8 

0.5763 

0.4237 

1.0 

0.6827 

0.3173 

1.2 

0.7699 

0.2301 

1.4 

0.8586 

0.1414 

1.6 

0.8904 

0.1096 

1.8 

0.9281 

0.0719 

2.0 

0.9545 

0.0455 

2.2 

0.9722 

0.0278 

2.4 

0.9836 

0.0164 

2.6 

0.9907 

0.0093 

2.8 

0.9949 

0.0051 

3.0 

0.9973 

0.0027 

3.2 

0.9986 

0.00137 

3.4 

0.9993 

0.00067 

3.6 

0.9997 

0.00032 

3.8 

0.9998 

0.00014 

4.0 

0.9999 

0.00006 

Prob[-Acr " 

x(t) ■ + Act] 

, PM 



—Xo 0 +\a x 
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x(f) 


° 'VV/V X V/'A^ 7 




the S(ojj against frequency curve. A typical spectral density function is shown in 
Figure 3.31. 

A random process whose spectral density is constant over a very wide frequency 
range is called white noise. If the spectral density of a process has a significant 
value over a narrower range of frequencies, but one which is nevertheless still wide 
compared with the center frequency of the band, it is termed a wide-band process 
(Figure 3.32). If the frequency range is narrow compared with the centre frequency 
it is termed a narrow-band process (Figure 3.33). Narrow-band processes frequently 
occur in engineering practice because real systems often respond strongly to specific 
exciting frequencies and thereby effectively act as a filter. 
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PART 5 

MECHANICS OF FLUIDS 

Peter Tucker 


3.10 INTRODUCTION 


Fluid is one of the two states in which matter can exist, the other being solid. In 
the fluid state the matter can flow; it will, in general, take the shape of its container. 
At rest a fluid is not able to sustain shear forces. 

Some “solids” may flow over a long period (glass window panes thicken at the 
base after a long time in a vertical position). The substances considered in fluid 
mechanics are those that are continuously fluid. 

Fluid mechanics is a study of the relationships between the effects of forces, 
energy and momentum occurring in and around a fluid system. The important prop¬ 
erties of a fluid in fluid mechanics terms are density, pressure, viscosity, surface 
tension, and, to some extent, temperature, all of which are intensive properties. 

Density is the mass per unit volume of the substance. Pressure is the force per 
unit area exerted by the fluid on its boundaries. Viscosity is a measure of the fluid’s 
resistance to flow and may be considered as internal friction. The higher the co¬ 
efficient of viscosity, the greater the resistance. Surface tension is a property related 
to intermolecular attraction in the free surface of a liquid resulting in the apparent 
presence of a very thin film on the surface. The meniscus at the intersection of a 
liquid and its container wall and capillarity are further examples of intermolecular 
attraction. 

Temperature is more relevant to thermodynamics than to fluid mechanics. It 
indicates the state of thermal equilibrium between two systems or, more loosely, 
the level of thermal energy in a system. 


3.11 FLUID STATICS 


Pressure at a Depth 

The variation of pressure p and depth h in a fluid of density p is given by 



(3.5) 


Most liquids are assumed to be of constant density p. In such a liquid the pressure 
at a depth h below a free surface is given by 

P = Po + Pgh (3.6) 

where p 0 is the pressure above the free surface. 

For gases equation (3.5) may be solved only if the relationship between p and 
h is known. A typical case is the atmosphere, where the relationship may be taken 
as polytropic or isothermal, depending on the altitude. Tables relating the properties 
of the atmosphere to altitude are readily available as the International Atmosphere. 

3.39 
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Pressure Measurement 

Pressure may be expressed as a pressure p in Pa, or as a pressure head h in m of 
the fluid concerned. For a fluid of density p, p = pgh. There are various instruments 
used to measure pressure. 

Manometers. Manometers are differential pressure-measuring devices, based on 
pressure due to columns of fluid. A typical U-tube manometer is shown in Figure 
3.34(a). The difference in pressure between vessel A containing a fluid of density 
p A and vessel B containing fluid of density p B is given by 

Pa~ Pb = Pb8 z b + (P m - P B )8 h ~ Pa8 z a ( 3 - 7 ) 

where h is the difference in the levels of the manometer fluid of density p m and 
Pm > Pa an d p,„ > p B . If p A — p B = p, then the difference in pressure head is 

Pa - Pb = Z B - Z A + (^ - 1 ) h (3.8) 

Pg \P ) 

If p m < p A and p m < p B then an inverted U-tube manometer is used as shown in 
Figure 3.34(b). In this case the pressure difference is 





FIGURE 3.34 



Downloaded from Digital Engineering Library @ McGraw-Hill (www.digitalengineeringlibrary.com) 
Copyright ©2004 The McGraw-Hill Companies. All rights reserved. 

Any use is subject to the Terms of Use as given at the website. 

























































MECHANICAL ENGINEERING PRINCIPLES 


MECHANICAL ENGINEERING PRINCIPLES 


3.41 


Pa ~ Pb ~ Pa8Z a + ( Pb - pjgh - p B gZ B 
and if p A = p B = p then the difference in pressure head is 


(3.9) 


Pa ~ Pb 
Pg 


~ Z A ~ Z B + h 



(3.10a) 


of if p m « p (the manometer fluid a gas and A and B containing liquid), 


Pa ~ Pb 

Pg 


~ Z A ~ Z B + h 


(3.10b) 


The accuracy of a U-tube manometer may be increased by sloping one of the 
legs to increase the movement of the fluid interface along the leg for a given 
difference in vertical height. This may be further enhanced by replacing the vertical 
leg by a reservoir and the inclined leg by a small-bore tube (Figure 3.34(d)). 

Another method is to increase the cross-sectional area of the ends of the legs 
(or one of the legs), as shown in Figure 3.34(c), so that a small movement of the 
free surfaces in the enlarged ends results in a large movement of the surface of 
separation. 

Dial Gauges. Most pressure dial guages make use of a Bourdon tube. This is a 
curved tube with an oval cross section. Increase in pressure causes the tube to 
straighten, decrease makes it bend. The movement of the free end turns a pointer 
over a scale, usually via a rack-and-pinion mechanism. The scale may be calibrated 
in the required pressure units. 

Diaphragm Gauges. In these gauges the pressure changes produce a movement 
in a diaphragm which may be detected by a displacement transducer, or by the 
output from strain gauges attached to the diaphragm surface. 

Piezoelectric Transducers. A piezoelectric crystal produces a voltage when de¬ 
formed by an external force. This induced charge is proportional to the impressed 
force and so the output can be used to supply a signal to a measuring device which 
may be calibrated in pressure units. 

Fortin Barometer. Barometers are used to measure the ambient or atmospheric 
pressure. In the Fortin barometer a column of mercury is supported by the atmos¬ 
pheric pressure acting on the surface of the mercury reservoir. The height h of the 
column above the reservoir surface, usually quoted as millimeters of mercury (mm 
Hg), may be converted to pressure units p 0 by 


Po ~ Pgh ~ 13-6 X 9.81/7 
= 133.42/7 Pa 


(3.11) 


Aneroid Barometer. In this device the atmospheric pressure tends to compress an 
evacuated bellows against the elasticity of the bellows. The movement of the free 
end of the bellows drives a pointer over a dial (or a pen over a drum graph) to 
indicate (or record) atmospheric pressure variations. 


Downloaded from Digital Engineering Library @ McGraw-Hill (www.digitalengineeringlibrary.com) 
Copyright ©2004 The McGraw-Hill Companies. All rights reserved. 

Any use is subject to the Terms of Use as given at the website. 




MECHANICAL ENGINEERING PRINCIPLES 


3.42 


SECTION THREE 


Buoyancy 

When a body is immersed in a fluid the difference in pressure over the depth of 
the body produces a displacement force on the body. The first recognition of this 
is attributed to Archimedes. 

Displacement Force. The buoyancy or displacement force F H on a body fully or 
partially immersed in a fluid is equal to the weight of the volume of the fluid 
equivalent to the immersed volume of the body (the weight of the displaced volume 
V D of the fluid): 


(3.12) 


F b = pgV D 


This buoyancy force acts vertically upwards through the centroid of the displaced 
volume, which is known as the center of buoyancy ( B ). If the buoyancy force is 
equal to the weight of the body then the body will float in the fluid. If the weight 
of the body is greater than the buoyancy force then the body will sink. If the 
buoyancy force is greater than the weight of the body then the body will rise. 

In a liquid, for example, a body will sink until the volume of liquid displaced 
has a weight which is equal to that of the body. If the body is more dense than the 
liquid then the body will not float at any depth in the liquid. A balloon will rise in 
air until the density of the air is such that the weight of the displaced volume of 
air is equal to the weight of the balloon. 

Stability of a Floating Body. Figure 3.35 shows bodies in various stages of equi¬ 
librium. A body is in stable equilibrium if a small displacement produces a restoring 
force or moment as for the ball in the saucer in Figure 3.35(a) or the floating bodies 
in (d) and (g). A body is in unstable equilibrium if a small displacement produces 
a disturbing force or moment as for the ball in Figure 3.35(b) or the floating bodies 
(e) and (h). A body is in neutral equilibrium if a small displacement produces no 
force or moment as for the ball in Figure 3.35(c) or the floating bodies in (f) and 


(i). 


For a partially immersed body, the point at which the line of action of the 
buoyancy force F B cuts the vertical center line of the floating body in the displaced 
position is known as the metacenter ( M ). For a floating body to be stable M must 
lie above the body’s center of gravity, G. If M lies below G the body is unstable; 
if M lies on G the body is in neutral equilibrium. The distance GM is known as 
the metacentric height. The distance of the metacenter above the center of buoyancy 



(3.13) 


where I second moment of area of the body at the water line (liquid surface) about 
its central axis normal to the direction of displacement. (See Table 3.5.) 

Period of Oscillation of a Stable Floating Body. A floating body oscillates with 
the periodic time T of a simple pendulum of length k 2 /GM, where k is the radius 
of gyration of the body about its axis of rotation. The periodic time is given by 



(3.14) 
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(a) 



V//5//////A 

(c) 



(d) (e) (f) 


Liquid 



(g) 


<h> (i) 


FIGURE 3.35 Stability. 


3.12 FLUID FLOW 


Definitions 


Continuity. For almost all analysis, a fluid is considered to be a continuum, that 
is, with nondiscontinuities or cavities in the flow stream. Cavitation, two-phase flow, 
“bubbly” flow, etc. are special cases with nonstandard relationships. 

Therefore for a continuum, by considering the flow through an elemental cuboid 
the continuity equation in three dimensions may be shown to be 


7 (PCt) + 7 (p<A) + 7 (pv z ) = 0 

dX dy dz 


(3.15) 

where v x is the fluid velocity in the x direction, etc. For a fluid of constant density 

(3.16) 


dv x dv dv. 

— + — + — = 0 
dX dy dz 


That is, the velocity of an incompressible fluid flow cannot increase in all three 
directions at the same time without producing discontinuity or cavitation. 
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TABLE 3.5 Second Moments of Area 



For two-dimensional flow: 


dv x dv 
— + — 
dx dy 


= 0 


(3.17) 


For one-dimensional flow the continuity equation may be linked with the conser¬ 
vation of mass, which states that for steady flow conditions mass flow rate, m, is 
constant throughout a flow system: 


m = pAv (3.18) 

where A is the cross-sectional area normal to the direction of flow. 


Circulation F. Circulation is defined as the line integral of the tangential velocity 
around a closed contour: 


r = 



(3.19) 


T is positive if the closed contour is on the left. 
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Vorticity £ Vorticity is defined as the circulation per unit area, and by considering 
the circulation around the element in Figure 3.36(a) it can be shown that 


T = dVy _dV, 
dxdy dx dy 


(3.20) 


Rotation w. Rotation is defined as the instantaneous mean angular velocity of 
two mutually perpendicular lines in a plane of the flow field. By considering the 
angular velocities of the two lines OA and OB in Figure 3.36(b) it can be shown 
that 


2 \dx dy ) 

or the rotation is equal to half the vorticity. 


(3.21) 


u+-l^ 

dy 2 



£ Ux 
dx 2 


(a) Vorticity 
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Stream Lines. The stream line is a line drawn in a flow stream which is every¬ 
where tangential to the direction of flow. A family of stream lines may be described 
mathematically by a stream function i ft, where iji fn(jc,v). Each stream line has the 
same function with a value of i// peculiar to that line. 

Stream Tubes. Since a line has no thickness, there can be no flow along a stream 
line. The stream tube is a concept introduced to enable flow along a stream line to 
be studied. It is a tube of infinitely small cross-section with a stream line as its 
axis. 


Energy. Energy is the stored form of heat and work. The basic concepts applied 
in fluid mechanics are: 


• The conservation of energy 

• That energy is transferred only as heat or work 

• That energy in a fluid flow system is stored only as internal energy, kinetic energy, 
or potential energy 


Other forms of energy (electrical, magnetic, chemical, etc.) may have to be taken 
into account in some circumstances, but are not usually included in general fluid 
mechanics relationships. 

Enthalpy and entropy need to be considered for gas flow analysis. The basic 
energy-flow equation is the steady-flow energy equation: 

Q + W = m\ (h + y + gZ j (3.22) 


where Q 
W 
h 

Z 

v 


the rate of heat transfer 

the rate of work transfer (power) 

the specific enthalpy (if e is the specific internal energy, p the pressure 
and p the fluid density, then h = e + (pip)) 
the height above some datum 
the mean velocity of flow 


Specific means “per unit mass.” For non-steady flow conditions, either quasi-steady 
techniques or the integration of infinitely small changes may be employed. 

Momentum. Momentum is the product of mass and velocity (mv). Newton’s 
laws of motion state that the force applied to a system may be equated to the rate 
of change of momentum of the system, in the direction of the force. The change 
in momentum may be related to time and/or displacement. In a steady flow situ¬ 
ation the change related to time is zero, so the change of momentum is usually 
taken to be the product of the mass flow rate and the change in velocity with 
displacement. Hence, the force applied across a system is 


F = thAv 


(3.23) 


where Av is the change in velocity in the direction of the force F. 

For flow in two or three dimensions the resultant force may be obtained by 
resolving the forces in the usual way. The flow round an expanding bend shown in 
Figure 3.37 is a typical example. The force in the x direction, F x , and the force in 
the y direction, F , are given by 
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FIGURE 3.37 Expanding bend. 


F x = p 1 A 1 + mv 1 - [p 2 A 2 + mv 2 )c os 0 (3.23a) 

F y = — (/V^2 + mv 2 )s in 0 (3.23b) 

from which the resultant force F R and its angle of inclination a can be found: 

= [FI + F]}^ (3.24) 

and 


a = tan 1 



(3.25) 


In pipe flow the pressure forces pA must not be forgotten. 


Equations of Motion 

Application of the momentum equation in three dimensions to an irrotational, in- 
viscid fluid flow leads to the Euler equation: 

Dv 1 

- 7^7 = - V P + sV/i (3.26) 

Ut p 

which for steady flow along a stream tube becomes: 

- dp + vdv + gdZ = 0 (3.27) 

P 

Integration for a constant-density fluid gives: 

2 

- + 7 - + gZ = constant (3.28) 

P 2 

These energy per unit mass terms may be converted to energy per unit weight 
terms, or heads, by dividing by g to give: 
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— + ^—I- Z = constant (3.29) 

Pg 2 g 

which is the Bernoulli (or constant head) equation. 

These equations are the generally more useful simplifications of the Navier- 
Stokes equation: 


Dv 

Dt 


P B 


Vp + V{w(Vu + Ve)} 


where B is the body force and e the rate of expansion. 


(3.30) 


Incompressible Pipe Flow 

Flow Regimes. The two major flow regimes are laminar and turbulent. Laminar 
flow may be fairly accurately modelled mathematically. The fluid moves in smooth 
layers and the velocity is everywhere tangential to the direction of motion. Any 
perturbations are quickly dampened out by the fluid viscosity. 

In turbulent flow the mathematical models usually need to be empirically mod¬ 
ified. Viscous damping may not be sufficient to control the perturbations, so that 
the fluid does not move in smooth layers and the instantaneous velocity may have 
components at an angle to the direction of motion. 

The ratio of inertia forces to viscous forces in a fluid flow is known as Reynolds ’ 
number (Re). In a pipe diameter D, with a fluid of density p and dynamic viscosity 
tj flowing with velocity v, Reynolds’ number Re = pDu/rj. 

A high value of Re > 2300 indicates relatively low damping, predicting turbulent 
flow. A low value of Re < 2000 indicates relatively high damping, predicting 
laminar flow. These values were suggested in an historical experiment by Osborne 
Reynolds. 

Pipe Losses (friction). Liquids (and gases under small pressure changes) flowing 
through pipes usually behave as incompressible fluids. Within the flow there is a 
relationship between the shear stress in the fluid and the gradient of the change of 
velocity across the flow. In most light liquids and gases, the relationship approxi¬ 
mates to the Newtonian one: 



(3.31) 


where ris the shear stress in the fluid, dvldy the gradient of the velocity distribution 
across the pipe and 17 the dynamic viscosity. 

The viscosity of the fluid produces not only the velocity variation across the 
flow but also a loss of useful energy along the pipe usually regarded as a friction 
loss. The force associated with this loss of energy appears as a shear force in the 
fluid at the pipe wall. A relationship between the shear stress at the pipe wall t 0 
and the friction coefficient, f is: 

\ PV 2 f (3.32) 

where v is the average flow velocity. 
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For use in pipe flow problems with viscous fluids the Bernoulli equation (3.29) 
may be adapted to include a head loss term, h L . Applied between two positions (1) 
and (2) in a pipe in a gravitational field of acceleration g, the head equation gives: 


Pi u? Pi vl 

- + ^ + Zi=- + ^ + Z 2 
Pg 2 g pg 2 g 


+ h. 


(3.33) 


where the head loss term h L is the loss of energy per unit weight of fluid flowing. 

Note that if a pump, say, is introduced between (1) and (2) an energy gain per 
unit weight term h w , equivalent to the output of the pump written as a head, should 
be added to the left-hand side of the equation to give 

— + + Z! + h w = ^ ^ + Z 2 + h L (3.34) 

Pg 2 g pg 2 g 


The relationship used to determine the head loss in a pipe depends on the flow 
regime in operation as well as the type and surface finish of the pipe wall. 

A mathematical analysis of laminar flow may be used to obtain an expression 
for the head loss along a pipe in terms of the fluid properties, pipe dimensions and 
flow velocity. Relating the pressure change along a length, L , of pipe of diameter, 
D, to the change in shear force across the flow produces Poiseuille’s equation: 


K = 32 


■qvL 

~PgD i 


(3.35) 


If the flow regime is turbulent, then the relationships in the flow cannot be easily 
described mathematically, but the head loss may be derived by equating the shear 
force at the pipe wall to the change in pressure force along the pipe. This gives 
the D’Arcy equation: 


K = 


4fLv 2 _ 
D 2 g 


(3.36) 


This relationship may also be established using dimensional analysis. 

Unfortunately, the friction coefficient, / is not a constant but depends on the 
type of flow and the roughness of the pipe walls. There are general relationships 
between f and Re which may be expressed as equations of varying complexity or 
as charts. For smooth pipes: 


1 

V? 


4 log 10 (2 R e Vf) 


1.6 


For rough pipes with a roughness size k this becomes: 


1 

Vf 



+ 3.48 


(3.37) 


(3.38) 


The Colebrook and White equation is a general or universal friction equation: 


= 3.47 - 4 log 10 

Vf 


(- + ^7=) 

\D ReVf / 


(3.39) 


It is, however, usually more useful to obtain values of / from a chart such as Figure 
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3.38. {Note: the value of / used in American equations for head losses is four times 
that used in the United Kingdom, so if values of f are obtained from American 
texts they should be moderated accordingly or the corresponding American equation 
used.) 

An empirical relationship widely used in water pipe work is the Hazen—Williams 
equation, usually written as: 

(h \ 0 ' 54 

v = 1.38 Cm 063 (3.40) 

where m is the ratio of the cross-sectional area of flow to the wetted perimeter 
known as the hydraulic mean diameter and C is a coefficient which depends on 
the condition of the pipe wall. 

Pipe Losses {Changes in Section ). When a fluid flows through a sharp (sudden) 
change in the cross-section of a pipe, energy is dissipated in the resulting turbulent 
eddies at the edge of the flow stream, producing a loss of head (or energy per unit 
weight). If the flow is from a smaller area to a larger one (sudden enlargement), 
the head loss is 


K 


(Vl ~ v 2 f 
2 g 


(3.41) 


When the flow is from a larger area to a smaller area (sudden contraction), the 
narrowed flow stream entering the smaller pipe is known as a vena contracta. The 
loss of head is assumed to be that due to a sudden enlargement from the vena 
contracta to the full area of the smaller pipe: 


K 


vj(± 
2 g \C c 



(3.42) 


The contraction coefficient, C c , is the ratio of the vena contracta area to that of the 
smaller pipe area. A typical value of C c is 0.6, which gives 


h 


L 



(3.43) 


which is also the head loss at the sharp entry to a pipe from a reservoir. Energy 
dissipation at changes in section, and pipe entry and exit, may be reduced by 
making the changes smooth and gradual, though this may be relatively costly. 

Other pipe fittings, such as valves, orifice plates and bends, produce varying 
values of head loss, usually quoted as a fraction of the velocity head {v 2 /2g). 


Pipe Networks. A system of pipes may be joined together either in series (one 
after the other) or parallel (all between the same point). The friction head loss 
across a system of pipes in series is the sum of the losses along each pipe individ¬ 
ually. The flow rate through each pipe will be the same. Using D’Arcy’s head loss 
equation: 


= ^W S + ^ 


LnV\ 

D 2 2g 


+ 4 f L " v " 

In D„2g 


(3.44) 


and 
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V = Ml = V 2 A 2 = • • • = v„A„ 


(3.45) 


If the system of pipes is connected in parallel the head loss across the system 
is equal to the head loss along any one of the pipes, when the flow has settled 
down to steady. The flow rate through the system is the sum of the flow rates along 
each pipe. Again using the D’Arcy equation: 




D, 2g 


D 2 2g 


L v 2 

= —~—— 

D„2g 


(3.46) 


V = Mi + v 2 A 2 + . . . uA (3.47) 

In addition, the rate of flow into each junction of a network, either in series or 
parallel, is equal to the rate of flow out of it. 

Pipe network problems are thus solved by setting up a number of such equations 
and solving them simultaneously. For a large number of pipes a computer program 
may be needed to handle the number of variables and equations. An example of a 
pipe network computer solution is given in Douglas et al. (1986). 


3.13 FLOW MEASUREMENT 


Pipe Flow 


Orifices and Nozzles (see Figure 3.39(a)). Another basic flow measurement tech¬ 
nique is to introduce some restriction into the flow passage and calibrate the re¬ 
sulting pressure changes against known flow rates. 

Often the restriction in a pipe is in the form of an orifice plate (a plate with a 
hole) or a nozzle. A simple application of the Bernoulli equation may be used for 
the design calculations, but it is always advisable to calibrate any measurement 
device in conditions as close to the required operating conditions as possible. 

Bernoulli and the continuity equations give the flow rate: 


V = 


C,,A 


o 


2-{p P - Po) ~r 
P [i - {AjA P n\ 


(3.48) 


where A a 

\ 

Pp 

Po 

Q 


the orifice (or nozzle throat) area 
the upstream pipe area 
the upstream pressure 

the pressure at the orifice or the nozzle throat 

a discharge coefficient which takes account of losses and contraction 
of the flow stream through the device 


Recommended orifice and nozzle dimensions, values of C d and methods of op¬ 
eration are contained in BS 1042. It is most important to place the orifice or nozzle 
so that its operation is not affected by perturbations in the upstream flow caused 
by valves, bends or other pipe fittings. 


Venturi Meters (see Figure 3.39(b)). The introduction of any restriction, particu¬ 
larly a sharp-edged orifice or nozzle, in a pipe will result in a loss of head. If it is 
required to keep this loss to a minimum, a Venturi meter may be used. The flow 
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(a) Orifice plate 





FIGURE 3.39 Flow meters. 


passage in a Venturi is gradually and smoothly reduced to a throat followed by a 
controlled expansion to full pipe section. In this way the head loss across the meter 
is greatly reduced, but the cost of producing a Venturi meter is much higher than 
that of an orifice. Equation (3.48) may be used to calculate the flow rate V, but the 
value of C d will now be approximately 0.98 for a well-designed Venturi meter. 
Again, BS 1042 should be consulted for recommended dimensions, values of C d 
and methods of operation. 


Rotameter or Gap Meter (see Figure 3.39(c)). If, somewhere within the system, 
it is acceptable to tolerate flow up a vertical section of piping, then a rotameter or 
gap meter may be used. This instrument depends on the balancing of the weight 
of a rotating float in a tapered glass tube with the drag forces in the annular passage 
surrounding the float. The drag forces depend on the flow rate and the correspond¬ 
ing area of the annulus. As the flow rate increases, the annulus area which will 
produce a drag force equal to the weight of the float also increases. Therefore the 
float moves up the tapered tube until the annulus area is such that the forces again 
balance. As the flow rate decreases the float descends to a reduced annulus area to 
again achieve a balance of forces. 
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SECTION THREE 


Velocity Meter. These are devices that measure velocity and not flow rate directly. 
Pitot and Pitot-static tubes are examples of such velocity-measuring instruments, 
making use of the pressure difference between the undisturbed flow stream and a 
point where the flow velocity is zero. They consist of two concentric tubes bent 
into an L shape as in Figure 3.39, with the outer tube joined to the inner at the toe 
of the L, at O. This end is usually spherical with a hole through to the inner tube. 
The undisturbed flow is assumed to be in the region of the holes round the periphery 
of the outer tube at X. The velocity is assumed zero at the spherical end presented 
to the flow, at O. 

The flow velocity, v may be calculated by applying Bernoulli’s equation between 
the two points O and X to give 



(3.49) 


P 


where p a is connected to O via the inner tube to the tapping at A, p x is connected 
to X via the outer tube to the tapping at B and C„ is a coefficient to cater for losses 
and disturbances not 

accounted for in Bernoulli’s equation. C v is often taken to be unity. The pressure 
difference may be measured using a manometer and then written into equation 
(3.49) as a head, h, to give 


0.5 



(3.50) 


As usual, it is advisable to calibrate the tube and obtain a calibration curve or an 
accurate value for C„. 

Care should be taken when a pitot-static tube is used to measure pipe flow, since 
the velocity will vary across the pipe. As a rough guide to the flow rate, the max¬ 
imum velocity, which is at the center of the pipe, may be taken to be twice the 
average velocity. Alternatively, the velocity at half the radius may be taken to be 
equal to the average velocity in the pipe. For an accurate evaluation the velocity 
distribution curve may be plotted and the flow rate through the pipe found by 
integration. This may be approximated to by dividing the cross-section into a series 
of concentric annuli of equal thickness, measuring the velocity at the middle of 
each annulus, multiplying by the corresponding annulus area, and adding to give 
the total flow rate. 

Current meters, torpedo-shaped devices with a propeller at the rear, may be 
inserted into pipes. The number of rotations of the propeller are counted electrically. 
This number together with coefficients peculiar to the propeller are used in empir¬ 
ical equations to determine the velocity. 

Velometers, varied anemometers, and hot wire anemometers are not usually used 
to measure the velocities of incompressible fluids in pipes, and will be discussed 
below under Gas Flow. 


3.14 BOUNDARY LAYER FLOW 


When a fluid flows over a solid boundary there is a region close to the boundary 
in which the fluid viscosity may be assumed to have an effect. Outside this region 
the fluid may be assumed inviscid. The viscous effect within the region is evidenced 
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by a reduction in velocity as the boundary is approached. Outside the region the 
velocity is constant. The region is known as a boundary layer. 

It is usual to assume that at the solid surface the fluid velocity is zero and at 
the boundary layer outer edge it is equal to the undisturbed flow velocity v s . This 
defines the boundary layer thickness S. (In practice, S may be taken to be the 
distance from the boundary surface at which the velocity is 99% of the undisturbed 
velocity, or 0.99 v s .) 

When a flow stream at a velocity v s passes over a flat plate, the boundary layer 
thickness S is found to increase with the distance x along the plate from the leading 
edge. Near the leading edge the flow inside the boundary layer may be assumed to 
be laminar, but as x increases the flow becomes turbulent and the rate of increase 
of S with x also increases, as shown in Figure 3.40. 

Within even a turbulent boundary layer there is a narrow region close to the 
plate surface where the flow is laminar. This is known as the viscous sublayer and 
has thickness S h . The reduction in velocity across the boundary layer is associated 
with a shear force at the plate surface, usually known as the drag force. 

Application of the momentum equation produces Von Karman ’s momentum in¬ 
tegral, in which the drag force per unit width, F D , becomes 



(3.51) 


where v is the velocity within the boundary layer at a distance y above the plate 
surface. (The integral 



may be defined as the momentum thickness ( 6 ) and the integral 



as the displacement thickness ( S *) so that 


f 



FIGURE 3.40 Boundary layer. 
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F d = pv 2 s e (3.52) 

In order to solve the Von Karman integral equation (3.51) or equation (3.52) it 
is necessary to know the value of 8 and the relationship between v and y, the 
velocity distribution. Both of these are dependent on each other and the flow regime, 
laminar or turbulent, within the boundary layer. 


Laminar Boundary Layers. A flat plate laminar boundary is normally assumed 
if Re x < 500 000. (Re x is Reynolds’ number based on x or pv s xl rj.) For laminar 
boundary layers various simplified velocity distribution relationships may be used, 
such as linear, sinusoidal, or cosinsusoidal. The generally accepted most accurate 
relationship is, however, that obtained by the reduction of a four-term polynomial, 
which gives 


£ = 3y _ 1 M 3 
v, 28 2 \8j 


(3.53) 


From this the shear stress at the plate surface, t 0 , may be found for Newtonian 
fluids: 


The shear force 




(3.54) 


(3.55) 


Substitution from equation (3.53) in equations (3.51) and (3.54) and equating 
F d from equations (3.51) and (3.55) leads via a separation of variables technique 
to 


^ = 4.64 Re x 05 (3.56) 

The drag force is usually quoted in terms of a drag coefficient, C D : 

2 F 

C D = ^ (3.57) 

pupc 

or 

F d = \ C dP v 2 x (3.58) 

By manipulation of the above equations an equation for C D for a laminar bound¬ 
ary over the whole length L of the plate: 

C D = 1 .29Re, 05 (3.59) 

is obtained. The drag force on the whole plate surface of area A is found from: 

F = \ C d v 2 A (3.60) 
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Turbulent Boundary Layers 


For Re x 500 000, a flat plate boundary layer is assumed to be turbulent. In a 
turbulent boundary layer the velocity distribution is often written in a power form: 


v 



(3.61) 


The index n varies between 6 and 9, depending on Re x . 

Because of the presence of the laminar sublayer, the turbulent regime is not 
continuous down to the plate surface, and (dv/dy) 0 does not give a useful result. 
The equation used for t 0 is 


T o = 0.0225 pv 2 R s 025 (3.62) 

based on work on smooth pipes by Blasius. 

Taking n — 7 and using the same techniques as for laminar boundary layers 
gives: 

- = 0.37 Re- 02 (3.63) 

and 

C D = 0.072 Re L 02 (3.64) 

This result assumes that the turbulent boundary layer obtains over the whole length 
of the plate to L. 

Prandtl suggested a more realistic expression which takes into account the pres¬ 
ence of a laminar boundary layer near the leading edge: 

C D = 0.074 Re L 02 - 1700 ReZ 1 (3.65) 

This may be used for 5 X 10 5 < Re, < 10 7 . For 10 7 < Re, < 10 9 Schlichting 
(1960) suggests a logarithmic velocity distribution and 

C D = O^log^e,)- 2 - 58 = 3.9l{lnRe L r 2 - 5S (3.66) 

Again, equation (3.60) may be applied to find the drag force on the whole plate. 


Viscous Sublayers 

The analyses above assume that the plate surface is smooth or at least hydraulically 
smooth. A surface is regarded as hydraulically smooth if the average roughness 
height k is less than the laminar sublayer thickness S b . For a turbulent layer with a 
velocity distribution power index of n = 1/7, the laminar sublayer thickness at a 
point at a distance x along the plate from the leading edge is given by 

^ = 199 ReZ 01 (3.67) 

o 

Thus S h may be compared with the roughness height, k, if the boundary layer 
thickness, 8, is known. 
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SECTION THREE 


3.15 PRESSURE TRANSIENTS 


Pressure transients may cause damage to piping systems such as aircraft fuel supply 
systems. It can be much more serious on a larger scale, where high-pressure rises 
over short periods may cause severe damage. Similar effects due to valve closure 
can be analysed on different levels of sophistication. 

The simplest is the rigid column theory, which assumes that the fluid is incom¬ 
pressible, and that the valve is closed relatively slowly. 


Slow Valve Closure 

When a fluid flowing through a pipe with a velocity v 0 undergoes a change in 
velocity there is an associated change in pressure. Equating the force due to the 
pressure change to the rate of change of momentum during closure gives the re¬ 
sulting pressure rise A p over a length of pipe L\ 



(3.68) 


The solution to this equation depends on a knowledge of the relationship between 
v and t (the valve closure rate in terms of the flow velocity). 

Equation (3.68) is only applicable to relatively slow valve closure rates in which 
the closure time should not be less than 2L/C (where C is the speed of sound in 
the fluid). 


Time to Establish Flow 

The rigid-column theory is also often used to calculate the time required to establish 
flow in a pipe on opening a valve. The theory implies that the time required to 
fully establish the flow is infinite and so the time t to achieve 99% of the final 
velocity v 0 is usually accepted: 



(3.69) 


where H is the supply head to the pipe entrance. The time t x required to reach x% 
of the final velocity is given by 



(3.70) 


Rapid Valve Closure 

When a fluid is brought to rest instantaneously from a velocity of v 0 by the closure 
of a valve at the exit of a pipe of diameter D , there will be a relatively high pressure 
rise at the valve. If the valve closure time is less than 2(L/c) then the resultant 
pressure rise is as if it were instantaneous; c is the speed at which the pressure 
wave travels through the fluid, which is the sonic velocity. 
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On such a rapid valve closure the kinetic energy of the flow is converted into 
strain energy in both the pipe material and the fluid (even liquids are acknowledged 
as compressible in this context). The resulting pressure wave is transmitted through 
the fluid away from the valve as shown in Figure 3.41. The pressure rise produced 
is 


A P = pcv 0 (3.71) 

For a fluid of bulk modulus G, in a pipe of wall thickness x, of a material with 
a Young’s modulus E and Poisson’s ratio cr, the velocity of the pressure wave is 


c = 



(3.72) 



if longitudinal stress is small compared to hoop stress. 


(3.73) 



Downloaded from Digital Engineering Library @ McGraw-Hill (www.digitalengineeringlibrary.com) 
Copyright ©2004 The McGraw-Hill Companies. All rights reserved. 

Any use is subject to the Terms of Use as given at the website. 



























MECHANICAL ENGINEERING PRINCIPLES 

3.60 SECTION THREE 


Equation (3.72) is often written as c = [G e l p\ 05 , where G e is the effective bulk 
modulus of the fluid and pipe combinations: 


G e 


I R. 

G + ~Ex 


(1.25 


o) 


(3.74) 


If the valve is at the entrance to the pipe, then rapid valve closure results in a 
rarefaction (pressure drop) at the valve. In other words, the pressure change is 


A P = ~pcv o 


(3.75) 


The Progress of a Pressure Wave 

Assuming no friction and no cavitation in the fluid, the progress of a pressure wave 
along a pipe between a valve and a reservoir following valve closure is as shown 
in Figure 3.41. The fluid in the pipe is successively brought to rest by the passage 
of the pressure wave. 

At a time t = Lie after valve closure, the pressure wave reaches the reservoir. 
The whole of the fluid in the pipe is at rest at a pressure p = p 0 + A p, which at 
the reservoir end instantaneously drops to reservoir pressure p 0 . The resulting pres¬ 
sure wave travels along the pipe towards the valve and the fluid at the higher 
pressure in the pipe flows towards the reservoir at its initial velocity v 0 . 

At t = 2L/c, the situation is the same as for a rapid closure of a valve down¬ 
stream of the flow, producing an instantaneous pressure drop to p 0 - A p and a 
rarefaction which travels towards the reservoir. The passage of the rarefaction suc¬ 
cessively brings the fluid to rest along the pipe. 

At t = 3 Lie, the rarefaction reaches the reservoir and the pressure instantane¬ 
ously rises to reservoir pressure p 0 . The resulting pressure wave travels towards the 
valve and fluid flows away from the reservoir at velocity v 0 . 

At t = 4L/c, the situation is the same as when the valve first closed at t = 0, 
and the cycle is repeated. 

In practice, friction quickly dampens out the pressure waves and cavitation re¬ 
duces the pressure decrease during the rarefactions. 

A typical plot of pressure against time at a valve following rapid valve closure 
is superimposed on the theoretical plot in Figure 3.42. 


3.16 GAS FLOW 


General Relationships 

The behavior of gases during processes involving thermal energy interactions and 
exchanges fits more properly into a study of thermodynamics. However, if only the 
flow mechanics are considered, the thermal and temperature effects may be re¬ 
stricted to those mainly relating to pressure and density. Applications of this ap¬ 
proach are found in aircraft and rocket engines. 

The most straightforward approach is to consider zero thermal energy transfer 
(heat transfer) to or from the fluid, or adiabatic flow. If, in addition, the changes in 
the fluid’s properties are assumed to be reversible, then the flow becomes isentropic 
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FIGURE 3.42 Pressure versus time at a valve. 


and the relationship between pressure and density the simple and well-known one 
for an ideal gas: 


— = constant 
P y 


(3.76) 


Also, the ideal gas law: 


= RT (3.77) 

P 

applies. Other useful relationships are: 

1. The ratio of specific heats y = c p /c v 

2. The gas constant R — c p — c u 

3. The universal gas constant R a — R X molecular mass 

= 8.3143 kJ kg m ', K > 

The terms stagnation or total temperature T 0 and pressure p 0 are often applied 
as the datum temperature and pressure of a fluid flow, even when stagnation con¬ 
ditions (zero velocity) do not exist in the particular situation under consideration. 
In gas flow the relationships between T 0 and the temperature T and p 0 and the 
pressure p at some point in the flow is often given in terms of the Mach number 
( M ), the ratio of the flow velocity v to that of sound c, i.e., 


(3.78) 

(3.79) 


M = 


c 


(3.80) 


and 
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c = 

In these terms T 0 and p 0 may be 

T 0 = 

Po = P 


SECTION THREE 


= [yRTf 


found from Euler’s equation to be 


1 + 




1 + 


(y - 1 )M 2 


yty-i) 


(3.81) 


(3.82) 

(3.83) 


Flow in Ducts and Nozzles 

Ducts. The analysis of gas flow in ducts is based on the Euler equation (3.26) 
and the one-dimensional continuity equation (3.18). Consideration of the differential 
forms of these equations will demonstrate that for subsonic flow (M < 1) the 
velocity will increase as the cross-sectional area of the duct decreases (in the con¬ 
verging entrance to a convergent/divergent nozzle, for example). For supersonic 
flow (M >1) the velocity will increase as the cross-sectional area increases (in the 
diffuser of the convergent /divergent nozzle). 

The properties of the fluid at a position in the flow stream where the local Mach 
number is unity are often denoted by a superscript * ( p *, p*, T*) and used as a 
datum, so that 


M* = 1; and v* = c* = [yi^T*] 0 - 5 


(3.84) 


The ratios of the properties at any position in the flow stream to those at the * 
position are: 


1 . 


P_ 

p* 


2 . 

3. 


P_ 

P* 

T_ 

p* 


4. 


5. 


A_ 

A* 


7+1 

_2 + (y - 1 )M 2 
y + 1 

2 + (y - 1) M 2 
y + 1 

2 + (y - 1 )M 2 


y/(y-l) 




M 

J_ 

M 


y + 1 

2 + (y - 1 )M 2 

2 + (y - 1 )M 2 
y + 1 


1/2 


(y+ l)/2(y—1) 


(3.85) 

(3.86) 

(3.87) 

(3.88) 

(3.89) 


For air the ratios may be calculated by substituting y = 1.4, or obtained from 
published tables and charts (Houghton and Brock 1961). 


Nozzles. A nozzle is an example of a duct with a smoothly decreasing cross- 
sectional area, followed in some cases by an increasing area (convergent/divergent 
nozzle) (see Figure 3.43). Since the velocity in the throat (minimum cross-section) 
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FIGURE 3.43 Convergent/divergent nozzle. 


3.63 


is often sonic, the approach velocity may be negligible. The throat velocity v t and 
the exit velocity v e are found by applying the Euler equation (3.26) between the 
upstream (entry conditions sufhx 0 ) and throat and exit, respectively: 


v, = 


2 Po 
T “ 1 Po 



(3.90) 


For v e the pressure term p e replaces p, in equation (3.90). The mass flow rate through 
the nozzle is usually found at the throat by substituting v, in the mass flow equation 
(3.19) to give 


m = C d A t 


2 Po 

y ~ i Po 



(3.91) 


where C d is a discharge coefficient which depends on the nozzle design. For a well- 
designed nozzle C d will be close to unity. 

The mass flow rate will be the same at the exit as at the throat. It may be 
calculated from the exit conditions by substituting A e and P e for A, and p t , respec¬ 
tively, in equation (3.91). 

Nozzles are usually designed for maximum mass flow rate. This will occur when 
the throat velocity is sonic ( v, = c). The pressure ratio which produces this situation 
is known as the critical pressure ratio, given by 



(3.92) 


For many light diatomic gases such as air, where y is approximately 1.4, {pt t )l 

(Po)crit = 0.528. 

The throat area will be that which gives the required mass flow rate through the 
throat at sonic velocity for critical pressure ratio. The exit area will be that which 
gives the calculated exit velocity for the given mass flow rate at the exit conditions. 
For convergent nozzles the throat also becomes the exit. 
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If nozzles, orifices, or Venturi meters are used to measure gas flow rates through 
a pipe then the approach velocity may be signficant and the mass flow rate given 
by 


m = C d A,p 0 


[ 2(po/p 0 )[p,/po) 2/ ' - (p f /p 0 )^ +1)/ '' ] 1/ 

l (y - 1)[1 - (pJp 0 ) 2,r (p,/Po) 2 J 


(3.93) 


Shock Waves 


Under normal design conditions the flow in the nozzle downstream of the throat 
will be supersonic. The velocity of the gas at exit will depend on the external 
pressure p b (back pressure). If the back pressure is greater than the theoretical exit 
pressure p e (p b > Pc) then shock waves will be set up in the nozzle. These are 
discontinuities similar to standing waves in open-channel flow. The shock waves 
set up in such a way are normal shock waves, normal to the direction of flow. 

If P b <Pe then the expansion will continue outside the nozzle (over-expansion). 

If conditions upstream of a normal shock wave are denoted by suffix 1 and 
downstream by suffix 2, then it can be shown that the product of the up- and 
downstream velocities is equal to the square of the sonic velocity at M — 1: 

v x v 2 = C* 2 (3.94) 


and since M l > 1, then M 2 < 1. Also, 
1. M 2 = 

Pi 

' Pi (7 +1) (y + 1) 

3 = b = (y + 1 W\ 

Pi v 2 2 + (y - l)Mj 


2 + (y - 1)M 2 
_2yM\ - (y - 1) 

2yM\ y — 1 


(3.95) 

(3.96) 

(3.97) 


These are known as the Rankine-Hugoniot relationships. Values for air may be 
obtained by putting y = 1.4 or by the use of published tables (Houghton and Brock 
1961). 

The strength of a shock wave may be defined as the ratio of the pressure change 
across the wave to the upstream pressure, or in terms of the upstream Mach number: 


Shock wave strength 


Pi ~ Pi 

Pi 



1.167(M 2 - 1), for air. 

(3.98) 


Oblique shock waves, at an angle f3 to the upstream flow direction, are produced 
when a supersonic gas flow is turned through an angle 6 by an obstruction such as 
an aircraft’s nose, wing or tail, the inside walls of a duct, etc. The relationships 
between the up- and downstream Mach numbers and the angles /3 and 0 are pub¬ 
lished in tables and charts (Houghton and Brock 1961). 

In some cases both the up- and downstreams will be supersonic and subsequent 
shock waves produced, for example, at the leading and trailing edges of a wing. 
The effects of such shock waves produced by aircraft in flight, say, may be noted 
at ground level as sonic booms. 
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Gas Flow Measurement 

Gas flow rates through ducts will normally be measured using devices and tech¬ 
niques similar to those used for incompressible fluids, namely orifice plates, venturi 
meters and nozzles. However, for gases the flow rate is usually quoted as a mass 
flow rate. Equations (3.90) and (3.92) may be used with orifices and Venturis as 
well as nozzles. Relevant values of C d for each device will be found in BS 1042, 
in addition to operational advice. 

When pitot-static tubes are used to measure gas flow velocities equation (3.49) 
may be acceptable for low flows with low pressure differences. At high velocities 
the compressibility must be taken into account and equation (3.90) used with stream 
conditions at X replacing those at the throat t. 

For accurate velocity measurement with little disturbance to the flow hot wire 
anemometers may be used. The resistance of an electrically heated wire is related 
to the temperature of the wire, which in turn is related to the velocity of the fluid 
flow past the wire. The wire resistance measured on a bridge may be calibrated 
against a known velocity, to give either direct readout or (more usually) a calibration 
curve. The fine wire of the anemometer is suceptible to fluid contamination. 

Other velometers and anemometers (depending on the relationship between the 
speed of rotation of a set of blades and the velocity (or speed) of the gas flow) 
may be used in very large cross-sectional ducts or to measure wind speed in the 
open air. They may depend on the rotation to generate a small electrical current, 
which can be calibrated as a speed, or the number of revolutions may be inserted 
into an empirical formula. A typical example is the three-vaned meteorological 
anemometer. Systems are now widely available for the measurement of velocity 
over a two dimensional field. Examples are laser doppler anemometry and particle 
image velocimetry, both of which require the seeding of the measured flow. 


3.17 IDEAL FLUID FLOW 


The concept of using idealized conditions to establish the shape of the mathematical 
models of real situations is common in engineering science studies. These models 
may then be modified to accommodate observed relationships, for application to 
real situations. 

Ideal fluid (or potential ) flow is such a concept. It may be used to set up flow 
patterns in the region of a flow stream outside the boundary layers. The combination 
of ideal flow and the boundary layer effects may be used to predict the performance 
of a real situation, so long as the limitations of both are recognized. The fluid is 
assumed to be inviscid and the flow steady, continuous, and irrotational. This means 
that there are no cavities or discontinuities in the flow stream, and that the fluid 
particles do not rotate about their own axes, even though the flow may be circular. 

The continuity equation (3.16) applies, and may be modified to 


^ + ^ = 0 
dx dy 


(3.99) 


if required for two-dimensional flow. For irrotation in two dimensions equation 
(3.21) becomes 


= 0 

dy dx 


(3.100) 
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The Stream Function. From the definitions of the stream line and stream function 
the equation to a stream line may be shown to be 


v x dy - v y dx = 0 (3.101) 

and since the stream function ip is an equation which describes a family of stream 
lines, then, for example. 


if/ = 2x — y 

represents a family of parallel straight lines with a variable intercept ip. 

For unit thickness in the z direction the volumetric flow rate V between two 
stream lines 1 and 2 is 


V = i//i - if/ 2 


or 


dV = dip 

Stream functions may be superposed so that if 

ip x = fn(jr,y) describes flow pattern A 


(3.102) 

(3.103) 


and 


ip 2 = fn(jr,y) describes flow pattern B 

then ip — ipi + >p 2 describes the flow pattern produced by the combination of A and 
B. 

The x and y components of the flow velocity v are given by 


v 


x 


dip 

—; and v 
dy 


dip 

dx 


(3.104) 


or 


dip dip 

v x = —; and v v =- 

* dy y dx 

depending on sign convention. In polar coordinates the components are 

— dip 


Radial velocity v r = 


rdd 


dip 

Tangential velocity v„ = — 


(3.105) 

(3.106) 


The Velocity Potential 

In a gravitational field there is a property the change in which is independent of 
the path of the change: potential energy. In a continuous, irrotational flow field 
there is also a property the change in which is independent of the path of the 
change. This property is the velocity potential (<p). It can be shown that 
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or 


d<j) 

; and v v 
dx 


d<f> 

dy 


Vg 


d(j) 

rdd’ 


d<f) 

and v, =- 

f dr 


(3.107) 


(3.108) 


Lines of constant (f> are known as velocity potential lines with an equation 

v x dx + v y dy = 0 (3.109) 

which intercept the stream lines at right angles to form an orthogonal network of 
characteristic pattern for each flow field. 

From equations (3.104) and (3.107) it can be seen that 


dl/r 

dy 


d<fi 

dx’ 


and 


dl/j 

dx 


d<f> 

dy 


(3.110) 


which are the Cauchy-Riemann equations. In such flow fields the Laplace equations 
for (f> and ijj must both be satisfied: 

V 2 0 = 0; and V 2 i//= 0 (3.111) 

Because the fluid is ideal, the Bernoulli equation (3.27) may be readily applied 
between points in the field, both along and across the stream lines. 


Flow Patterns 

Examples of simple flow patterns are shown in Figure 3.44 with the equations to 
their stream functions and velocity potentials. A source is mathematically a point 
at which fluid appears and flows radially outwards. A sink is a negative source at 
which fluid flows radially inwards to disappear at a point (similar to the plug hole 
in a domestic sink, where, however, there is a vortex superposed to produce a spiral 
vortex or whirlpool). A vortex is flow in concentric circles with no radial flow. A 
doublet is the superposition of a source and a sink of equal strength m, initially a 
distance 2 a apart brought infinitely close together so that the product of their 
strength and the distance between them remains a constant k. k is the strength of 
the doublet and is equal to 2 am. 

This is a mathematical concept which is apparently impractical but yields a 
useful flow pattern. It is often used in combination with other simple patterns. 


Modeling 

Since there can be no flow across a stream line and the fluid is assumed inviscid, 
any stream line in a flow pattern may be replaced by a solid surface with no effect 
on the rest of the pattern. A stream line forming a closed contour may be replaced 
by a solid body to model the flow pattern around a body of the same shape. This 
provides a method of writing mathematical models to describe the flow streams 
around various shapes. 

The limitations of the model must be understood, and the effects of rotation and 
viscosity particularly considered, when applying the analysis to real situations. For 
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FIGURE 3.44 Simple flow patterns. 
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example, the flow round a cylinder may be modeled by a combination of doublet 
and parallel flow. The drag and lift forces calculated by integrating the resulting 
pressure forces at the cylinder surface appear to be zero. This is obviously not 
correct. 

In the real situation the viscosity of the fluid produces a boundary layer at the 
cylinder surface, which, because of curvature, separates from the surface to form a 
wake. The presence of the wake disturbs the downstream flow pattern and the form 
drag force is a consequence. Viscosity also introduces a surface friction drag. 

The shedding of vortices from the cylinder surface into the wake also produces 
alternate positive and negative lift forces, which are not predicted by the ideal flow 
analysis, although telegraph wires may often be observed vibrating in the wind. 


3.18 CONCLUSION 


Fluid mechanics is often regarded as an empirical subject which makes use of 
formulae based only on observed experimental results. This misconception is further 
compounded by the extensive use of coefficients (discharge especially) to account 
for effects which are difficult to model mathematically. However, almost all fluid 
mechanics equations in common use are based on the conservation of energy, the 
conservation of momentum or the fact that the rate of change of momentum may 
be equated to an applied force, usually a pressure force. The experimental checking 
and empirical amendment to derived formulae is just good engineering practice. 

This section does not provide rigorous derivations of the various equations 
quoted. Some appreciation of such derivations may be required in order to establish 
the limitations and modifications necessary for the application of the equations, 
especially to nonstandard situations. For this, standard textbooks on fluid mechanics 
as listed below should be consulted. 
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PRINCIPLES OF 
THERMODYNAMICS 

Dennis H. Bacon 


3.19 INTRODUCTION 


Thermodynamics is concerned with energy transfers in processes. Two modes of 
transfer are recognized: work (transfer) and heat (transfer). Before proceeding, the 
terminology used in thermodynamics must be defined. 

The particular part of the working substance under consideration is called the 
system, and this is separated from the surroundings by a boundary. In the closed 
or nonflow system the mass of working substance is constant, but in the open or 
flow system there is a mass flow rate across the boundary. Some processes in 
reciprocating plant may be considered by nonflow analysis but in steam plant, for 
example, most are considered by flow analysis. The state of a system is defined by 
the properties (pressure, temperature, etc.). Properties are normally expressed spe¬ 
cifically (i.e., per unit mass) to enable charts or tables to be used. The state of 
simple substances can be described by two independent properties, but complex 
ones such as mixtures need more definition. A change of state is achieved by a 
process which is idealized as reversible with no losses. Reversible processes can 
be described by mathematical equations and enable analysis to be made to give 
answers for ideal situations. Real processes have losses and are described as irre¬ 
versible, and the ideal results are multiplied by a coefficient or efficiency (based 
on measurement or experience) to predict real performance. 


3.20 THE LAWS OF THERMODYNAMICS 


The First Law of Thermodynamics 

This is a law of energy conservation. When applied to a process we write 
Q — W = A E or q — w — Ae 

where Q is the heat transfer (kj) or q is the specific heat transfer (kj/kg), W is the 
work transfer (kj) or w is the specific work transfer (kj/kg), and A E is the energy 
change (kj) or Ae is the specific energy change (kj/kg). 

The change symbol A means final value minus initial value. A E embraces all 
forms of energy, but in the nonflow process it is usual to find that the only signif¬ 
icant change is in the internal energy ( U,u ) and we write the nonflow energy equa¬ 
tion 


Q — W = AU or q — w = An 
For the steady flow system we write 
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{ V 2 \ ( V 2 

Q ~ W x = ;«A I h + — + gz I or q - w x = A I h + — + gz 

where Q and W x are the energy transfer rates and m is the steady mass flow rate 
across the boundary (in and out), A h is the change in specific enthalpy ( h = u + 
pv), W 2 I2 is the change in specific kinetic energy, and A gz is the change in specific 
potential energy. The suffix x is used on the work transfer to denote that this is the 
useful work from the system as the flow work is included in the enthalpy term. In 
flow problems it will also be necessary to use the continuity equation 

m = pAV 

where p is the density and A is the area normal to the velocity V. Analysis of 
nonsteady flow may also be made, in which case energy terms to allow for the 
storage of energy in the system will be added. 

Warning: A sign convention for work and heat is built into the equations above. 
Positive work means work obtained from the system and positive heat means heat 
put into the system. Care should be taken to be clear about the symbol V, which 
may appear as velocity or volume in many equations. 

In order to allow continuous energy transfers a cycle is defined in which a series 
of processes brings the working substance back to the initial state so that the cycle 
can be repeated continuously. If we apply the first law to a cycle it follows that A E 
is zero and 


2 Q = 2 w 

cycle cycle 


The Second Law of Thermodynamics 

It might be thought that the first law of thermodynamics permits all the heat transfer 
to a cycle to be returned as work transfer, but unfortunately the second law places 
restraints on the achievement of this desirable situation. The restraint takes the 
practical form of demanding that some of the heat transfer to the cycle must be 
rejected as a heat transfer to a lower temperature. Thus when we build a heat engine 
it has to exchange heat with (at least) two reservoirs in order to produce work 
(Figure 3.45). Since work is the objective, the amount produced per unit heat input 
is vital information and we define the thermal efficiency of a heat engine as 

Net work transfer from the cycle W 
'hihermai | [cat transfer to the cycle Q x 

Since the first law states Q 1 — Q 2 = W we see that efficiency is less than unity. 

The second law makes further investigations and determines the maximum pos¬ 
sible efficiency of a heat engine using reversible isothermal processes to transfer 
heat from two reservoirs as 


Thermal maximum ' * min' 'max) 

This efficiency is known as the Carnot efficiency and is not attainable due to losses. 
It is also found that constant temperature processes, except during phase change, 
are not practical and real processes of heat transfer take place at approximately 


Downloaded from Digital Engineering Library @ McGraw-Hill (www.digitalengineeringlibrary.com) 
Copyright ©2004 The McGraw-Hill Companies. All rights reserved. 

Any use is subject to the Terms of Use as given at the website. 



MECHANICAL ENGINEERING PRINCIPLES 

3.72 SECTION THREE 



FIGURE 3.45 A heat engine. 


constant volume or constant pressure. One positive product of the second law is 
that it tells the engineer that thermal efficiency will be increased by making the 
maximum cycle temperature as high as possible (a materials constraint) and by 
making the minimum cycle temperature as low as possible (ambient conditions). 

The losses in a real cycle are due to internal fluid friction and the necessity of 
having a temperature difference to cause a heat transfer. The larger the temperature 
difference, the greater the losses. The fluid friction losses in a work-producing 
process are defined by the process efficiency 

Actual work produced 
^/process id ea l work produced 

which is inverted for work-absorbing processes (compression). The most common 
application of this efficiency is in steady flow adiabatic processes. Ideally, these 
are processes with no heat transfer which are often used as models for real processes 
in which the heat transfers are negligible compared with the work transfers (turbine 
expansion). In such processes the efficiency is known as the isentropic efficiency 
because the ideal adiabatic process has constant entropy. 

Losses due to fluid friction and losses due to heat transfer across finite temper¬ 
ature differences are found to result in an increase in the value of the entropy that 
would be expected in a reversible process. Thus, an expected increase would be 
larger and an expected decrease would be smaller. It is not easy to define entropy 
except mathematically. In practical use as the abscissa of charts it enables work 
transfers in ideal adiabatic processes to be represented as vertical lines if enthalpy 
is used as an ordinate, and in this guise is a valuable visual method of presentation. 


3.21 THERMOECONOMICS 


When a more detailed study, see Kotas (1985) and Bejan (1982), of a flow process 
is made by the second law of thermodynamics it is found that specific entropy (J / 
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kgK) appears as part of a property known as availability. In a flow process we 
write b = h — TgS, in which b is the specific availability function, h is the specific 
enthalpy, T 0 is the temperature (absolute) of the surroundings, and .v is the specific 
entropy. The second law shows that the maximum work potential or exergy of any 
state in surroundings at state 0 is given by b — b 0 . Thus for a change of state in a 
flow process from 1 to 2 the maximum specific work obtainable is given by the 
exergy change, w Vmax = (b x — b 0 ) — (b 2 — b Q ) = (b 1 — b 2 ) = —A b. If we measure 
or predict by analysis the actual work achieved it is possible to determine numer¬ 
ically the lost work or irreversibility in the process. If engineering plant is to be 
designed to the best advantage it is clear that processes should be chosen to min¬ 
imize this loss. The lost work may be associated with costs and we move into the 
developing field of thermoeconomics. Clearly, this is a complex subject but it is 
important in that it unites thermodynamics with costs and can help in the design 
of long-life expensive plant, such as aircraft engines. 


3.22 WORK, HEAT, PROPERTY VALUES, 
PROCESS LAWS AND COMBUSTION 


To deploy the laws of thermodynamics outlined above we need more information. 
To perform simple cycle analysis the data below are vital. 


Work 

In a nonflow process work transfer can be determined from w = fpdv. The math¬ 
ematical relation for the process is known as the process law ( qv ). In most flow 
processes used in engineering cycles the adiabatic approximation is used so that 
the steady flow energy equation, neglecting changes in kinetic and potential energy, 
gives 


w x = A/r 


Heat 

This is usually an unknown quantity and is found by the application of the energy 
equation. As stated earlier, many processes are approximately adiabatic so that heat 
transfer is zero and in others heat transfer is obtained from combustion data or, if 
a heat exchange process, by heat exchanger efficiency. 


Property Values 

These are found in tables or from charts for common substances. Computer for¬ 
mulations are now widely available. 


Process Laws 

This is a particularly important step in thermodynamic analysis because an idealized 
reversible process has to be chosen to represent as closely as possible the real 
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process in order to calculate energy changes. When the working substance is a gas, 
it is convenient in an elementary analysis to use perfect gas laws with the process 
calculation. These are 


pv = RT or pV = mRT 
u — c v \T and h = c p \T 

where c„ and c p are the specific heat capacities at constant volume and constant 
pressure, respectively, which are related as follows: 

c p ~ c„ = R and c„/c„ = y 

where R is the specific gas constant and y is the isentropic index. 

Ideal processes commonly used are constant pressure, constant volume, constant 
temperature (which for a perfect gas becomes pv = constant) together with two 
other more general relations: the adiabatic process, pv k = constant (which for a 
perfect gas becomes pv y = constant), and the polytropic process, pv" = constant. 
The last process is a general relation between pressure and volume which is used 
if none of the other clearly special cases are considered valid. Usually 1 < n < 
1.4. 

It is possible (by using the gas laws) in adiabatic and polytropic gas processes 
to rearrange the relations to involve pressure and temperature or temperature and 
volume to yield very useful relations: 

1± = (p± V^ 1)/T . h = h = (p±Y"~ 1Vn . T± = feV 1 
t 2 W ’ t 2 \vj ' t 2 \pj ’ t 2 \vj 

Processes may be represented on property diagrams to enable cycle visualization 
(Figures 3.46-3.51). .. 


Combustion 

To avoid involving complex chemical equations, engineers often use the calorific 
value of a fuel coupled with a combustion efficiency to estimate the energy transfers 
in combustion processes. Thus, the rate of energy input by combustion is 

E = rh f - CV ■ Tj comb 

where m, is the fuel mass flow rate, CV the calorific value of the fuel, and i 7 comb 
the combustion efficiency. 


3.23 CYCLE ANALYSIS 


One example will be given of the simple analysis of the ideal Joule cycle for a gas 
turbine plant (Figure 3.51). The cycle consists of four flow processes described in 
Table 3.6 and analyzed by the steady flow energy equation. 

From the data in the table it can be seen that the specific work w = c p (T 2 - 
T 4 ) — c p (T 2 — TJ and the thermal efficiency 
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V 


FIGURE 3.46 Substance phases and definitions. 


h 


P 



s 

FIGURE 3.47 Steam processes on an h-s dia¬ 
gram. 
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FIGURE 3.48 Steam processes on a T-s di¬ 
agram. 



FIGURE 3.49 Gas processes on a p—v dia¬ 
gram. 


= c p (t 3 - t a ) - c p (t 2 - r t ) 

thermal ^ _ yy 

If allowance is made for the isentropic efficiency of the compression and expansion 
processes the cycle diagram is changed to show the associated entropy increases 
but the expressions for work and efficiency above are still valid with the changed 
values of T 2 and T 4 (Figure 3.52). These values are determined from the use of the 
reversible adiabatic process relation and the isentropic efficiency as 

T 2 - r, = - l)/ Vc 


and 
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5 

FIGURE 3.50 Gas processes on a T—s diagram. 



FIGURE 3.51 The Joule cycle. 


T, ~ T 4 = V ,T 3 (] - 1 /r<r u *) 


where r p is the cycle pressure ratio. 


Vc = 


7V 

T 2 


~ T\ 

~ T, 


and 


V, = 



17 , and 17 , being the isentropic efficiencies of compression and expansion. If these 
values are substituted into the work and thermal efficiency expressions they become 
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TABLE 3.6 


Process 

Description 

q — w x = h 

1 to 2 

Reversible adiabatic compression 
pv 7 = constant 

^T 

I 

E-h 

ii 

i 

-sT 

II 

i 

2 to 3 

Reversible constant pressure heat 
transfer to the cycle 

q = h 3 - h 2 = c[p(T 3 - T 2 ) 

3 to 4 

Reversible adiabatic expansion 
pv y = constant 

-w x = h 4 - h 3 w x = c p (T 3 - T 4 ) 

4 to 1 

Reversible constant pressure heat 
transfer from the cycle 

q = h 1 - h 4 = c p ( 7’, - T 4 ) 



FIGURE 3.52 The effect of isentropic process 
efficiency on the Joule cycle. 



FIGURE 3.53 The effect of pressure ratio in 
a cycle with fixed T m . lx and 7" lnlI1 . 
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Design 


range 


T3 

C 

CO 

JL 

cr 



FIGURE 3.54 Thermal efficiency and specific 
work transfer variation in a Joule cycle with al¬ 
lowance for isentropic process efficiency. 


more useful in that they involve the thermodynamically significant maximum and 
minimum cycle temperatures which are fixed by material and ambient conditions 
respectively, so that the only variable is the cycle pressure ratio (Figure 3.53). If 
the expressions are differentiated with respect to this pressure ratio it is possible to 
find the pressure ratio for maximum work and that for maximum efficiency. The 
cycle designer then has a choice, depending on the proposed application and Figure 
3.54 shows that it would be expected that the chosen ratio would fall between these 
two maxima. Obviously, this simple approach is not the complete answer to gas 
turbine cycle analysis, but it illustrates the use of the laws of thermodynamics, and 
similar work may be done for other plant cycles. 
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3.24 INTRODUCTION 


Whenever a temperature difference occurs there is an energy flow from the higher 
temperature to the lower. A study of heat transfer is concerned with the determi¬ 
nation of the instananeous rates of energy flow in all situations. We determine heat 
transfer rates in watts. These rates will be constant in situations where the temper¬ 
ature difference remains constant but variable (transient) when the temperature dif¬ 
ference varies either due to the heat transfers or to other energy changes such as 
internal chemical reaction. 

There are three modes of heat transfer: 

1. Conduction, which is of greatest interest in solid bodies but also occurs in fluids, 
where it is often overshadowed by convection 

2. Convection, which occurs in fluids when energy is transferred due to the motion 
of the fluid 

3. Radiation, which occurs between two systems at different temperatures which 
need not be in contact provided any intervening medium is transparent to the 
radiation. 

In practice, all three modes may occur simultaneously and it is necessary to 
draw up a balance at a boundary. For example, energy may be conducted to the 
surface of an electric storage heater and is then convected and radiated to the 
surroundings. Thus calculations can become complex, and in this particular case 
where energy is added at certain times this is a continuously varying situation. 
Three approaches to heat transfer will be discussed below: 

1. A simple method suitable for many estimations 

2. A more detailed appraisal of the field 

3. Comments on the use of computers 

3.25 BASIC PRINCIPLES OF HEAT TRANSFER 
(White 1984; Ozisik 1985; Kreith and Bohn 1986 ) 


Conduction 


Fourier’s law for conduction states: 


<7 = 



The thermal conductivity t(Wm 'K ') is a property of the material which varies 
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with temperature but for small temperature ranges is usually considered constant. 
Typical values are shown in Table 3.7. With constant k Fourier’s equation can be 
integrated for four common situations. 


Plane Surfaces. Integration gives 


Q 


k l2 A 




12 


(see Figure 3.55). The quantity (A x/kA) is known as the thermal resistance in 
KW Thermal resistances can be added in a similar way to electrical resistances 
so that for a multilayer plane surface there are a number of resistances in series 
(Figure 3.56). Thus we can write 


Q = - T 4)/S ^ 


Cylindrical Surfaces. For tubes it is more convenient to evaluate heat transfer 
rates per unit length, and integration gives 


Q = (T, 


T 2 )/ 


In r 2 /r l 
2t rk l2 L 


A 77 


In r 2 lr x 
2irk l2 L 


(see Figure 3.57), and in this case the thermal resistance is 


TABLE 3.7 


Substance 

Thermal conductivity, Wm 1 K ', at 20°C 

Aluminum 

204 

Iron 

52 

Water 

0.597 

Air 

0.026 (100 kPa) 

Glass wool 

0.04 



FIGURE 3.55 One-dimensional conduction 
through a single-layer plane wall. 


Multilayer plane surface 



FIGURE 3.56 One-dimensional conduction 
through a multilayer plane. 
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FIGURE 3.57 Conduction through a single- layer cylindrical surface; [/-value for a plane 
layer cylindrical surface. surface. 


/ in r 2 /r l \ 
\ 2Trk n L / 


in KW- 1 


For a multilayer tube, thermal resistances are added to give 


Q' = (F, 


T 4 )l E 


/ In (fou,er/finner) \ 

\ 27 TkL ) 


(see Figure 3.58). 



Fluid 

2 


Bulk 


temperature 

7f 2 


(b) 


FIGURE 3.58b {/-value for a cylindrical surface. 
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Convection 

The fundamental equation for convective heat transfer at a solid-fluid interface is 


Q" = he 


where 9 is the temperature difference between surface and fluid. The surface heat 
transfer coefficient /i(Wnf 2 K _1 ) is not a property of the fluid or the surface but 
depends on the flow rate, the fluid properties, and the surface shape. The coefficient 
has to be determined for each situation and can vary considerably (Table 3.8). 
Although the determination of h is crucial to convection calculations, it is an ex¬ 
tremely difficult process, and accurate prediction of convective heat transfer is not 
always possible. 

If we express the convection equation in a thermal resistance form suitable for 
plane surfaces, 


it can be seen that the thermal resistance is (2/nA)KW 1 . For tubular surfaces it is 
again more convenient to work per unit length, so that 



and the thermal resistance is (l/2i7r/t)mKW 1 . 


Overall Heat Transfer Coefficients 

A common heat transfer situation is a solid wall separating two fluids, and for this 
problem the thermal resistances for conduction and convection can be added to 
enable the heat transfer rate to be determined in terms of the two fluid temperatures. 
For a plane surface (Figure 3.58(b), 



For a tubular surface (Figure 3.59), 


TABLE 3.8 Range of values of surface heat 
transfer coefficient (W m 2 K ') 


Free convection Gases 

Liquids 

Forced convection Gases 

Liquids 


10 to 700 
100 to 10 100 


0.5 to 500 
50 to 2000 
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Fluid 1, Tf, 

Boundary layer, h<\ 

Double-layer tube 
Boundary layer, h 7 
Fluid 2, 7" f2 


1 _ 1 2 In outer/ r inner) 1 

6/' 2^/7, 2irk 2rrr 2 *2 


FIGURE 3.59 1/'-value for a cylindrical surface. 



O' = ('/V, - rj, 


1 


2irr l h 1 




ln(r, 


/ r. 

outer' ' inner / 1 


2 irk 


2tt r 2 h 2 


It can be seen that the added resistances may be inverted to give an overall con¬ 
ductance which is known as a {/-value or overall heat transfer coefficient. 

For a plane surface, 

I - 1 y — 1 

U h l + ^ k + h 2 

For a tubular surface, 

J_ = 1 + y / ln(rou, er /r in ner) \ 1 

U' 2m\h l \ 2-7 rk ) 2irr 2 h 2 

The heat transfer rate is then simply written 
For a plane surface, 

Q = UA(T fl - T h ) where A is the area 


For a tubular surface 


Q = U'l(T h - T fl ) 


where / is the length. 

The situation in which this technique is commonly used is in heat exchanger 
design. It should suffice for simple calculations provided suitable values of the 
surface heat transfer coefficients for convection can be obtained. (CIBS Guide 
1980). If the temperature difference is not constant then a mean value should be 
used. A suitable equation for a mean can be found in the heat exchanger section 
which follows. 
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Radiation 

Radiation is of central importance in space application where all rejected heat is 
radiated to space. The rate of energy emitted by an ideal black body is given by 
the Stefan-Boltzmann law, in which the absolute temperature (Celsius + 273) is 
raised to the fourth power: 


El = oT 4 

where cr is the Stefan-Boltzmann constant 5.67 X 10 s Wm 2 K 4 and the subscript 
b refers to the ideal black body. Real bodies emit less radiation, and the mono¬ 
chromatic emissivity is defined by 


F" 

■ C 'A 


The value of s varies with A and T because real bodies are selective emitters, but 
for simple calculations it is often assumed that emissivity is constant. The calcu¬ 
lations associated with this assumption are based on grey body theory, for which 
the rate of energy emission is given by 

E" g = eoT 4 

It would be unwise to estimate unknown emissivities, and measurements would 
need to be made unless suitable data could be found. 

Radiation incident on a body may be absorbed, reflected or transmitted. Thus 
we write a + p + r = 1 where a, p and r are the absorptivity, reflectivity, and 
transmissivity, respectively. Ideal black bodies absorb all incident radiation but real 
bodies do not. Gases are often assumed to transmit all radiation, but this is not 
always true, particularly with hydrocarbon combustion products and atmospheric 
transmission. Solids have a transmissivity of zero. With these simple ideas it is 
necessary to know the values of only a and p. It can be shown that a gray body 
has absorbtivity equal to emissivity, a = s. Thus, provided the transmissivity is 
zero, a knowledge of the grey body emissivity enables reflectivity to be determined, 
since 


p = 1 = £ 

The only simple radiation problem that can be solved with the simple approach 
above is that of a gray body in large surroundings (see Figure 3.60). The word 
‘large’ implies that radiation not incident on the grey body which will be incident 
on the surroundings and will therefore be reflected will not be re-incident on the 
gray body. Thus, the surroundings are effectively black. (This might be true with 
a linear size factor greater than 10.) In this simple case it can be shown that the 
heat transfer rate is 


q = e<tA(t 4 - r 4 ) 

where s is the emissivity of the body, A is the area of the body, T is the temperature 
of the body in K and T s is the temperature of the large (black) surroundings in K. 
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Surroundings 


FIGURE 3.60 A small gray body in large 
(black) surroundings. 


Simple Transient Problems 

If a body is being cooled or heated by convection or radiation and the thermal 
conductivity is large so that the rapid heat transfer rates within the body enable it 
to be assumed that the body temperature distribution is uniform, then the situation 
is known as a lumped capacity system. For such a system the complex methods of 
transient heat transfer are not required and a simple energy balance equation may 
be drawn up and integrated. The most common case is quenching, a convective 
boundary problem for which in time dt a small heat transfer SQ occurs when the 
body temperature changes from T by an amount dT (Figure 1.61). Thus, 


SQ = pc„VdT = -hA(T - T f ) 


where T, is the fluid temperature, A is the body surface area, V is the body volume, 
p is the body density, c p is the body specific heat, and h is the surface heat transfer 
coefficient. Integration gives 


- = e ~(hAr/pCpV) 

00 


where 6 0 is the initial temperature difference between fluid and body and 0 is the 
temperature difference at any future time t. The quantity ( pc p V/hA) may be re¬ 
garded as the time constant of the system. 



Fluid 


Body 


coefficient, h 

0 = (T — r,) 


Specific heat, C p 
Temperature, T 


FIGURE 3.61 A lumped-capacity system with 
convection. 


Downloaded from Digital Engineering Library @ McGraw-Hill (www.digitalengineeringlibrary.com) 
Copyright ©2004 The McGraw-Hill Companies. All rights reserved. 

Any use is subject to the Terms of Use as given at the website. 



MECHANICAL ENGINEERING PRINCIPLES 

MECHANICAL ENGINEERING PRINCIPLES 3.87 

3.26 ANALYSIS OF HEAT TRANSFER 


Conduction 

By considering the thermal equilibrium of a small, three-dimensional element of 
solid, isotropic material it can be shown that for a rectangular coordinate system 


dT 

dt 


d 2 T d 2 T d 2 TO'" 
to 2 dy- to J pc p 


where dT/dt is the rate of change of temperature with time, a is the thermal dif- 
fusivity of the material a = k/pc p and Q" is the internal heat generation rate per 
unit volume, which may be due, for example, to electric current flow for which 
Q'" = Ft, where i is the current density and r the resistivity. The solution to this 
equation is not easy, and numerical approximation methods are often used. One 
such method is the finite difference technique, in which continuously varying tem¬ 
peratures are assumed to change in finite steps. Consider the three planes shown a 
distance Ax apart (Figure 3.62). At 


and at 


so that at 


dT _ 7', - T 2 
to Ax 

dT T, ~ T , 

B — = — k -- 

to Ax 


d 2 T T + T 3 -2T 2 
dx 2 Ax 2 


Similarly, dT/dt may be written 



FIGURE 3.62 One-dimensional finite difference formulation. 
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- T 


A t 


where T n l is the temperature at layer n at time 1 and 7'„ 0 is the temperature at 
layer n at time 0. 

For steady state situations dT/dt = 0 and a two-dimensional plane surface will 
be used for illustration requiring a solution of 


d 2 T d 2 T~\ Q'" 
dx- dy 2 \ pc p 


Consider the surface to be divided by a grid (Figure 3.63). It is then found that 
the solution for any point in the plane for steady state conduction without heat 
generation is 

T, + T 2 + r 3 + T 4 - 4 T 0 = 0 

or with heat generation is 

a 2 Q”’ 

T, + T 2 + r 3 + T 4 - 4 T 0 =- f- 


At the boundary of the plane conditions are usually isothermal, in which case T = 
constant, or convective when an energy balance yields for a straight boundary (Fig¬ 
ure 3.64) 


— + T 
2 + 12 






= o 


where h is the surface heat transfer coefficient, T f is the fluid temperature, and a is 
the grid size. Similar expressions can be derived for corners, curves, etc. at the 
boundary. 

For a large number of grid points, a large number of simultaneous equations are 
obtained which can be solved by iteration or Gaussian elimination. Computer pro¬ 
grams may be used to advantage. The solution obtained will be the temperature 



FIGURE 3.63 The two-dimensional grid con- FIGURE 3.64 Convective surface nomencla- 
cept. ture. 
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distribution in the plane, and the heat transfer rates may be found at the boundary 
(Figure 3.65): 

With an isothermal boundary Q = "Zk{T m — r wall ) 

With a convective boundary Q = 'Zha(T f - T m ) 

For transient heating or cooling for which dT/dt =A 0, a one-dimensional illus¬ 
tration is used. The equation to be solved is 

dT _ (FT 
dt 01 dx 2 

when there is no heat generation (Figure 3.66).expressed in finite difference form 
this becomes 


Tn.i F T n+10 + T n0 2^ 

where F is the nondimensional grid size Fourier number F = a\t/a 2 . The only 
unknown in this equation is T n l , the temperature at layer n after one time interval 
At. Thus from a knowledge of the initial conditions successive temperatures in each 
layer can be found directly for each time interval. This is the explicit method and 
is used for tabular or graphical (Schmidt method) solutions. If F > 0.5 the solution 


7"mi Tmj 


Fluid 

Vs%//A / //A/s Wal1 

Body 


T, = constant 


Fluid 

T w aii = constant 

7fr77fi777%7s wall 


Body 


Tm-\ Tmj T m3 


(a) Convective boundary (b) Isothermal boundary 

FIGURE 3.65 Calculation of heat transfer rate. 




FIGURE 3.66 One-dimensional transient conduction formula¬ 
tion. 
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is unstable and in three dimensions the criterion becomes severe. The boundary 
conditions may be isothermal or convective and in the latter case the solution is 


T n ,i = 2 F • 7',, , ,, + r„ 0 [l - IF - IF ■ B] + 2F ■ BT { 


where B is the nondimensional grid Biot number, B = ha Ik. For this case the 
solution is unstable if (F + FB) > 0.5. The solutions obtained give the temperature 
distribution in the one-dimensional plane and the heat transfer is found at the bound¬ 
ary 



or from the temperature profile 


Q 2 m C p fl final ^initial) 


layer 


The stability problems of the explicit method can be overcome by the use of 
implicit methods for which there is no direct solution, but a set of simultaneous 
equations are obtained which may be solved by Gaussian elimination. A computer 
program may be used to advantage. A satisfactory implicit method is that due to 
Crank and Nicolson. The importance of a stable solution is that if the choice of F 
is limited then the grid size and time interval cannot be freely selected, leading to 
excessive calculations for solution. The implicit method releases this constraint but 
care is still needed to ensure accuracy. 

Although the finite difference method has been chosen for demonstration be¬ 
cause the method is easy to understand, most modem computer programs are based 
on the finite element technique. However, the mathematical principles are involved, 
and would not lend themselves to simple programming. Before the availability of 
computer software analytical solutions were obtained and presented as graphs of 
transient solutions for slabs, cylinders, and spheres. These graphs enable solutions 
for other shapes to be obtained by superposition methods. Such methods should be 
used to avoid or validate computer solutions. 

Warning: If fiber-reinforced materials are used in which the lay-up is arranged 
to give directional structural strength it will be found that the thermal conductivity 
has directional variation and the methods above will need considerable amendment. 


Convection 

A knowledge of the surface heat transfer coefficient li is essential in determining 
heat transfer rates. Fluid flow over a solid surface is a boundary layer problem, and 
the heat transfer depends on boundary layer analysis. This analysis may be by 
differential or integral approach, but solution is difficult and the modeling of tur¬ 
bulence is complex. Computer solutions based on numerical approximations may 
be used to advantage, but simple approaches have been used for many years and 
are still extremely useful. These methods are based on Reynolds’ analogy (modified 
by later workers) and dimensional analysis backed by experimentation. 

Convection may be free or forced. In forced convection it is found that the heat 
transfer coefficient can be included in a nondimensional relation of the form 


Downloaded from Digital Engineering Library @ McGraw-Hill (www.digitalengineeringlibrary.com) 
Copyright ©2004 The McGraw-Hill Companies. All rights reserved. 

Any use is subject to the Terms of Use as given at the website. 


MECHANICAL ENGINEERING PRINCIPLES 

MECHANICAL ENGINEERING PRINCIPLES 3.91 


Nu = 4>(Re,Pr) = constant • Re" ■ Pr b 

where Nu is the Nusselt number (Nu = hi Ik), Re is the Reynolds number (Re = 
pVl! pi), and Pr is the Prandtl number ( Pr = pr p /k). In these relations / is a rep¬ 
resentative length dimension (diameter for a pipe and some chosen length for a 
plate), V is the bulk or free stream velocity outside the boundary layer. The values 
of the constants a and h depend on whether the flow is laminar or turbulent and 
on the geometry of the situation, and are usually found by experiment. 

The determination of whether flow is laminar or turbulent is by the value of the 
Reynolds number; 

For plates. Re < 500 000, flow is laminar: Re > 500 000, flow is turbulent 

For tubes, Re < 2000, flow is laminar: Re > 4000, flow is turbulent 

(between these two values there is a transition zone). There are many relations to 
be found in texts which allow for entry length problems, boundary conditions, etc. 
and it is not feasible to list them all here. Two relations are given below which 
give average values of Nusselt number over a finite length of plate or tube in forced, 
turbulent flow with Mach number less than 0.3 using total plate length and diameter 
for representative length dimension. Care must be taken in any empirical relation 
to use it as the author intended. 

Plate: Nu = 0.036 \Re°*Pr 033 

In this relation fluid properties should be evaluated at the film temperature, T am = 
(^wall "P 7bulk)/2- 

Tube: Nu = 0.023Re 08 Pr 04 

In this relation fluid properties should be evaluated at the bulk temperature, 0.6 < 
Pr < 160 and (l/d) > 60. 

It should be noted that the index of Reynolds number of 0.8 is characteristic of 
turbulent flow; in laminar flow 0.5 is found. 

It must be emphasized that reference to other texts in all but these simple cases 
is essential to estimate heat transfer coefficients. It should also be pointed out that 
the values obtained from such relations could give errors of 25%, and a search of 
the literature might reveal equations more suited to a particular situation. However, 
an estimate within 25% is better than no knowledge, and is a suitable starting point 
which may be modified in the light of experience. 

For complex heat exchange surfaces such as turbine blade cooling in an aircraft 
engine, empirical information is usually presented graphically (on these graphs the 
nondimensional group St (Stanton number) may appear: 



RePr pVc p 


In free convection the relationship used is Nu — 4>(Pr ■ Gr), where Gr is the Grashof 
number, p 2 f}gdl 3 /p? in which /3 is the coefficient of cubical expansion of the fluid 
and 0 is a temperature difference (usually surface to free stream temperature). The 
transition from laminar to turbulent flow is determined by the product (Pr ■ Gr) 
known as the Rayleigh number, Ra. As a simple example, for plane or cylindrical 
vertical surfaces, it is found that 
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For Ra < 10 9 , flow is laminar and Nu = 0.59 (Pr ■ Gr) 025 

For Ra > 10 9 , flow is turbulent and Nu = 0.13(Pr • Gr) 1 ' 3 

The representative length dimension is height and the resulting heat transfer 
coefficients are average values for the whole height. Film temperature is used for 
fluid properties. Warnings similar to those given for forced convection apply to the 
use of these equations. 

Phase change convection heat transfer (condensing and evaporation) shows co¬ 
efficients that are, in general, higher than those found in single-phase flow. They 
are not discussed here but information may be found in standard texts. 


Radiation 

In space technology applications heat transfer systems are designed for conduction 
and radiation transfer only. When an emitting body is not surrounded by the re¬ 
ceiver, the spatial distribution of energy from the radiating point needs to be known. 
To determine this distribution the intensity of radiation i^ is defined in any direction 
4> as 



where dco is a small solid angle subtended at the radiating point by the area inter¬ 
cepting the radiation, dco = dA/r 2 (Figure 3.67) (the solid angle represented by a 
sphere is 4tt steradians). Lambert’s law of diffuse radiation states that i. = i n cos 
c t> where i n is the normal intensity of radiation which can be determined for black 
and gray bodies; 




FIGURE 3.67 Spatial distribution of radia- FIGURE 3.68 Heat transfer by radiation be- 
tion. tween two arbitrarily disposed gray surfaces. 
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Black i n crT 4 1 tt Gray i n scrT 4 / ir 

With this knowledge of the radiation intensity in any direction it is only nec¬ 
essary to determine the amount that any body can see of any other body to calculate 
the heat transfer rate. For this purely geometric problem mathematical analysis 
(Figure 3.68) suggests a quantity variously known as the geometric, configuration 
or shape factor, which is defined as the fraction of the energy emitted per unit time 
by one surface that is intercepted by another surface. The geometric factor is given 
by 


_ 1 f f cos</>! cos </> 2 dA x dA 2 
12 Aj Jai Ja 2 ttx 2 

It can be seen that A l F l2 = A 2 F 2l , a useful reciprocal relation. It is also clear 
that the equation will require skill to solve in some situations, and to overcome this 
problem geometric factors are available for many situations in tables or on graphs 
(Hottel charts). The charts can give more information than anticipated by the use 
of shape factor algebra, which enables factors to be found by addition, subtraction, 
etc. (Figure 3.69). 

Having established the intensity of radiation and the geometric factor, problems 
may be solved by an electrical analogy using the radiosity of a surface. Radiosity 
is defined as the total emitted energy from a gray surface: 

j" = E" s + pG" 

where j" is the radiosity and pG" the reflected portion of the incident radiation 
G". Since E" — eE" the net rate of radiation leaving a gray surface of area A 
becomes 


El - J" 

p/As 

which may be envisaged as a potential difference, El — j", divided by a resistance, 
p/AE. A similar geometric resistance of 1/AFcan be established to enable complete 
circuits to be drawn up. Thus for a three-body problem we may sketch the analo¬ 
gous electrical circuit (Figure 3.70) and apply Kirchhoff’s electric current law to 
each J" node to obtain three simultaneous equations of the form 

I ” _ j" in _ in in _ in 

In J 1 _j_ J 2 J 1 J 3 J \ _ q 

p 1 /A 1 E 1 1 /A | F| 2 1/A^ 

If there are more than three bodies sketching becomes complex and the equation 
above can be rearranged and generalized. For N surfaces ( / = 1 to A) there will 
be N equations, the rth of which (i = 1 to N) will be 

J'!-(I- £,) i F i jj; = e ,£" 

When j = i, F u will be zero unless the the surface is concave and can see itself. 
This set of N simultaneous equations may be solved by Gaussian elimination for 
which a computer program may be used. The output of the solution will be N 


Downloaded from Digital Engineering Library @ McGraw-Hill (www.digitalengineeringlibrary.com) 
Copyright ©2004 The McGraw-Hill Companies. All rights reserved. 

Any use is subject to the Terms of Use as given at the website. 







MECHANICAL ENGINEERING PRINCIPLES 



3.94 

Downloaded from Digital Engineering Library @ McGraw-Hill (www.digitalengineeringlibrary.com) 
Copyright ©2004 The McGraw-Hill Companies. All rights reserved. 

Any use is subject to the Terms of Use as given at the website. 








































MECHANICAL ENGINEERING PRINCIPLES 

MECHANICAL ENGINEERING PRINCIPLES 3.95 



values of j" and any required heat flux can be found. In the three-body problem, 
for example, the heat transfer from 1 to 2 


Q 12 


j" ~ J 2 

l/A t F 12 


or the total heat transfer from body 1 is 


p l /A 1 s 1 


Special cases are: 

1. If IV — 1 bodies are in large surroundings then for this Mh body j" = E" h . 

2. An insulated or refractory surface has no black body potential but contributes 
to the heat transfer by taking up an equilibrium temperature T given by crT 4 = 

j,, 

3. Radiation shield problems will show six resistances in series rather than the 
series-parallel circuits used previously. 


In all the discussion on radiation above no account has been taken of the selec¬ 
tive emitter for which emissivity is not constant. Additional techniques are required 
to solve these real problems and care should be exercised if widely varying emis¬ 
sivity is encountered. In new situations it will be necessary to determine emissivity 
by experiment before proceeding. No account has been taken of intervening media 
for which transmissivity is not unity. Gas absorption and radiation needs further 
information. For gases such as oxygen, nitrogen and hydrogen with symmetric, 
diatomic molecules, the above work is adequate, but asymmetric molecular struc- 
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tures cause problems. In particular, hydrocarbon fuel combustion products (H 2 0, 
C0 2 , CO, S0 2 ) are important in engineering calculations and account must be taken 
of their radiation properties. Solar radiation problems also need special considera¬ 
tion. 


Finned Surfaces 

In many heat exchange problems involving the determination of a 17-value it is 
found that the surface heat transfer coefficient on one side of the solid interface is 
much smaller than that on the other. The smaller coefficient will dictate the heat 
transfer rate achieved, and in order to overcome this problem fins may be added to 
this poor convection surface to increase the area for heat transfer. This problem will 
occur in liquid to gas exchangers on the gas side. The addition of fins will alter 
the flow pattern so that a new coefficient should if possible be determined. It is 
also possible that there may be variation of coefficient over the fin surface. 

Simple fin theory in which conduction along the fin is balanced with convection 
from the surface can be used to determine the temperature distribution and heat 
transfer rate of the fin. For example, when a long fin of constant cross-sectional 
area is examined (Figure 3.71) it is found that the temperature distribution is 

0 cosh(/n(/ — x )) 

6 0 cosh (ml) 


and the heat transfer rate is 


Q = mkA6 0 tanhfm/) 

where 8 is the temperature difference between fin and fluid and 9 0 is the difference 
at the fin root, m = hp/kA, p is the fin perimeter and A is the cross-sectional area. 

A fin efficiency is then defined to compensate for the varying temperature dif¬ 
ference along the fin as 


Vfin = 

which for the 


Actual heat transfer rate 

Heat transfer rate if whole fin were at the wall temperature 
simple case above is Tj fln = tanh (ml)/ml. Fins are usually fitted in 



So » T„ - Ti 

FIGURE 3.71 Simple rectangular 
fin nomenclature. 
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arrays and the efficiency of a fin system can be established in the form of an area 
weighted fin efficiency, 17 ' 


17 ' = %m/3 +1-/3 


where 


total fin area A fin 

total area including fins A 

The U value based on the enhanced area A is then 

1 1 _ v Ax 1 

U A ~ (A 1 /A)h 1 ~ ^ R + v 7;, 

For complete surfaces /3 is supplied by the manufacturer. 

When fins of more complex shape are used (tapered fins or annular fins) the 
cross-sectional area is not constant and fin efficiency data are obtained from graphs. 
Care should be taken in the interpretation of such graphs since the equations above 
may not agree with the definitions used for the graphs. 


Heat Exchangers 

It is possible to design a heat exchanger with the information above and obtain a 
basic idea of size and configuration of simple tubular structures. It is first necessary 
to realize that temperature differences change along a heat exchanger, and that flow 
may be parallel or counter in direction. The latter is to be preferred, as it leads to 
smaller sizes. To allow for the changing temperature difference a log mean tem¬ 
perature difference is used (Figure 3.72): 


ln( 6 t / 6 2 ) 


and to allow for varying flow patterns (which are neither counter nor parallel flow) 
graphs are available to give a factor F to modify the (usually) counterflow value 
of mean temperature difference. Thus, for any heat exchanger, 


If an estimated mean U -value for the surface is then determined the heat exchange 
equation is 

Q — UAF6 LMrD 

Thus the area is determined. There are many solutions to this equation to satisfy 
all the constraints which will include: 


Energy for each stream Q = ( mAh) hol = (mAfl) cold (A h = enthalpy change) 
Continuity of each stream til = pAV 
Heat transfer area A = mil 


The above equations will involve options with varying tube numbers and di- 
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iParsIlel or counter 

FIGURE 3.72 The mean temperature difference in heat exchang¬ 
ers. 


ameters which will affect the determination of heat transfer coefficients (and hence 
the (/-value), and a number of solutions will be obtained. The optimum choice will 
involve allowable pressure drops, velocities and exit temperatures. It is not a dif¬ 
ficult calculation, but, because of the choice, a computer program may be used. 
There are short cuts to this approach based on the interrelation between pressure 
drop and heat transfer (modified Reynolds’ analogy) to determine heat transfer 
coefficients and a method based on graphs of effectiveness ( E ), capacity ratio (c), 
and number of transfer units ( NTU ) is sometimes used. 

By whatever method a design is achieved, it will, unless the application is very 
simple, be necessary to consult a professional heat exchanger designer with expe¬ 
rience and full computer programs. 


3.27 USE OF COMPUTERS 


Computers may be used as an aid to heat transfer calculations at various levels of 
skill. At the simplest level, computer programs may be written to determine, for 
example, heat losses in buildings, heat inputs from pipes and radiators, etc. This 
will save repetitive calculations and build a small library of useful programs. To 
aid those whose heat transfer and computing skills are slight, Bacon (1989) gives 
a simple approach to problems. For more detail on the mathematics of finite dif¬ 
ference techniques in heat transfer Myers (1971) is useful. 

An alternative approach to numerical approximations in computer work is to use 
finite element methods. For heat transfer applications Myers (1971) is again useful. 
Very few people will find the need for finite element programming skills, as there 
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is an ever-increasing range of software aimed at the solution of many engineering 
problems, including heat transfer. Many of these are PC packages in which there 
is integration between design, drawing, manufacture and analysis of stress, vibra¬ 
tion, heat transfer, etc. 

Whenever numerical approximation techniques are used, large-scale analysis 
leads to considerable demand on data storage and computing time. It is therefore 
essential to do as much as possible with simple methods before becoming com¬ 
mitted to large-scale finite element packages. If it is decided that the use of such a 
package is necessary it is vital to be sure of the requirements of the problem, for 
volume of input data and output results will be large and unpalatable. To assist 
with this problem, preprocessing packages are used for data input and mesh gen¬ 
eration with graphic display, and post-processing packages are used for graphic 
display of the results. For example, the temperature distribution in a combustion 
chamber or piston displayed by color graphics enables easier identification of prob¬ 
lem areas. 

Finally, it must be emphasized that computing is not a substitute for understand¬ 
ing heat transfer problems. It is an aid to enable a more detailed investigation to 
be achieved and presented in a fashion to enable engineers to improve their designs. 


3.28 HEAT TRANSFER: NOMENCLATURE 

A area 

B Biot number 

c p Specific heat capacity at constant pressure 
E" Black-body emissive power 
E' g Gray-body emissive power 
F Fourier number, geometric factor, mean temperature factor 
g gravitational constant 
Gr Grashof number 
h surface heat transfer coefficient 
i. intensity of radiation in direction <fi 
i h normal intensity of radiation 
j” radiosity 
k thermal conductivity 
/ length 

m mass, fin parameter 
tit mass flow rate 
Nu Nusselt number 
p perimeter 
Pr Prandtl number 
Q heat transfer 
Q heat transfer rate 
Q' heat transfer rate per unit length 
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Q" heat transfer rate per unit area 
Q'" heat transfer rate per unit volume 
r radius 

Ra Rayleigh number 
Re Reynolds number 
St Stanton number 
T temperature 
t time 

U overall heat transfer coefficient per unit area 
U' overall heat transfer coefficient per unit length 
V velocity 

x, y, z rectangular coordinates 

Greek Letters 

a thermal diffusivity, absorbtivity 
;8 coefficient of cubical expansion 
A change in 
17 efficiency 
e emissivity 
A wavelength 
I x viscosity 
<f> angle 

p density, reflectivity 
a Stefan-Boltzmann constant 
S summation 
0 temperature difference 
r transmissivity 
<p solid angle 
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